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The environs of supermassive black holes are among the universe's most extreme phenom- 
ena. Understanding the physical processes occurring in the vicinity of black holes may 
provide the key to answer a number of fundamental astrophysical questions including 
^ the detectability of strong gravity effects, the formation and propagation of relativistic 

I ^^ I jets, the origin of the highest energy gamma-rays and cosmic-rays, and the nature and 

evolution of the central engine in Active Galactic Nuclei ( AGN) . As a step towards this 

^vj direction, this paper reviews some of the progress achieved in the field based on obser- 

►^ vations in the very high energy domain. It particularly focuses on non-thermal particle 

/"^v acceleration and emission processes that may occur in the rotating magnetospheres orig- 

I inating from accreting, supermassive black hole systems. Topics covered include direct 

I electric field acceleration in the black hole's magnetosphere, ultra-high energy cosmic 

^vj ray production, Blandford-Znajek mechanism, centrifugal acceleration and magnetic re- 

• connection, along with the relevant efficiency constraints imposed by interactions with 

C^*" matter, radiation and fields. By way of application, a detailed discussion of well-known 

O sources (Sgr A*; Cen A; M87; NGC1399) is presented. 



Keywords: Black Hole, Magnetosphere, Particle Acceleration, Radiation Mechanism, 
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1. INTRODUCTION 

C^ This review focuses on non-thermal particle acceleration and very high energy 

( VHE) emission processes that may occur in the vicinity of magnetized supermassive 
black holes. As such, active galaxies and, to some extent, extinct or dormant quasars 
are placed in the center of its attention. Although rotating, supermassive black holes 
have often been considered as putative sources for the energization of ultra-high en- 
ergy (UHE) cosmic rays, non-thermal magnetospheric models have recently gained 
an additional impetus with the detection of very high energy gamma-rays from 
the (non-blazar) radio galaxy M87, located ^ 16 Mpc away in the Virgo cluster of 
galaxies EEEl While former observations of rapid VHE flux variations (on timescales 
of ~ [1 — 2] days) already indicated a small size of the 7-ray emitting region, the 
location of this region remained ambiguous based on the VHE results alone. Ad- 
ditional high-resolution (^ 50 Schwarzschild radii) radio VLBA imaging has now 
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shown that the gamma-ray flaring activity is accompanied by an increase in the 
radio flux close to core (see Fig. IT]), indicating that the required energetic charged 
particles may in fact be accelerated in the very vicinity of the central black hole™ 
At present, further evidence is needed to strengthen this inference. Nevertheless, 
magnetospheric models where the relevant non-thermal processes are considered to 
occur at the base of a rotating black-hole-jet magnetosphere (e.g., refsEEEl), have 
emerged as attractive candidates. If verified by further observations, VHE gamma- 
ray observations could provide a fundamental diagnostic of the most violent region 
in Active Galactic Nuclei (AGN). 

This review is structured as follows: Sect. 2 gives a phenomenologically-orientated 
introduction into the close environment of supermassive black holes. In Sect. 3, 
particle acceleration scenarios that draw on rotating, large-scale poloidal magnetic 
fields (i.e., the magnetosphere) are presented. Sect. 4 discusses possible interactions 
of accelerated particles with ambient matter, radiation and fields, that are expected 
to limit achievable particle energies. Applications to concrete sources as potential 
VHE gamma-ray and UHE cosmic ray emitters are discussed in Sect. 5. 



2. SUPERMASSIVE BLACK HOLES AND MAGNETIC FIELDS 

2.1. Evidence and Mass Range 

Most, if not all galaxies, are nowadays believed to harbor a supermassive black hole 
(BH) with 10^ to lO^-^A/0 at their centeJUl (for reviews see also refsHliiJ). Initially 
hypothesized in order to account for the huge power output seen from extragalactic 
radio sources such as quasars,Q2l the presence of supermassive black holes has now 
directly been probed by a variety of means, e.g., through the velocity dispersion of 
stars, water maser line emission or the kinematics of nearby ionized gas. The most 
prominent example is probably the center of our own galaxy where near-infrared 
imaging revealed proper motion of stars which increases with a Kepler law down to 
separations of less than five light days from the compact radio source Sgr A*, thus 
providing strong evidence for the presence of a black hole of mass ~ 3 x IO^MqISES 
Several important correlations for the black hole masses in galaxies have been es- 
tablished within recent years, including (i) a relation between the black hole mass 
and the luminosity/mass of its host galaxy bulge, e.g., Mbh — 2 x lO^'^A/buigc, 
and (ii) a relation between the black hole mass and the host galaxy stellar velocity 
dispersion a, log(MBH/Mo) = (8.21 ± 0.06) + (3.83 ± 0.21) log(cr/200 km/s) (with 
intrinsic scatter of 0.22 ± 0.06 dex).'i3For the nearby (z=0.034) gamma-ray blazar 
Mkn 501 (HBL source), for example, a — 291 ± 13 km/s has been found, ^- sug- 
gesting a central black hole mass of ~ 3 x lO^M©. Similarly, observations of the 
giant elliptical radio galaxy NGC 1399 in the centre of the Fornax cluster [d ~ 20 
Mpc) suggest a ~ 320 km/s,'^^ indicating the presence of a black hole with mass 
exceeding ^ 5 x IQ^Mq (see also reLl^SJ). 
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Fig. 1. Combined VHE gamma-ray, X-ray and radio light curves for the radio galaxy M87 in 
2007-2008. A VHE gamma-ray flare (Fig. A) is followed by an increase of the radio flux (Fig. C) 
close to the black hole. Knot HST-1, located in the jet ^ 100 pc away, is in a low X-ray state 
(Fig. B) and therefore unlikely to be the source of the observed VHE gamma-rays. Details are 
as follows: (A) VHE -ray flux data {E > 0.35 TeV, nightly average) including H.E.S.S., MAGIC, 
and VERITAS. The inlay shows a zoomed version of the rapid flaring activity (with timescales as 
low as ~ [1 — 2] days) in February 2008. (B) X-ray (Chandra) flux for the nucleus and the knot 
HST-1. (C) Radio (43 GHz VLBA) flux data. The shaded horizontal area indicates the range of 
radio fluxes from the nucleus before the 2008 flare. The radio flux of the outer jet regions does not 
change substantially; most of the observed increase results from the region around the nucleus. 
Figure adapted from Acciari et al.l^ 



2.2. Rotating Black Holes 

Supermassive black holes residing in the centers of galaxies may be driven into 
rotation by prolonged accretion of angular momentum or as a result of merger 
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Fig. 2. Schematic drawing of the structure of a rotating Kerr black hole. 

events. 



2.2.1. Characteristic properties 

The angular momentum J for a (rotating) Kerr black hole of mass M is usually 
expressed in terms of the dimensionless spin parameter a < 1, J — a Jmax, where 
the maximum value Jmax is given by 

J - ^^ (1) 

'-'max — ■ \^ } 

C 

SO that an extreme (maximally spinning) Kerr black hole is characterized by a spin 
parameter 

' ' ^1. (2) 



GA'P/c J„, 

For a non-rotating (Schwarzschild) black hole we have J = a = 0, and the event 
horizon scale rn coincides with the Schwarzschild radius 

2GM ,^ f M \ 

^s = ^-3xl0"(^^cm. (3) 

For a rotating Kerr black hole the (spherical) event horizon surface is located at 

(l + Vl-a^) , (4) 






with < a < 1, and where r^ — rs/2 is often called the gravitational radius. Hence, 
for a Schwarzschild black hole r}j — r^ and for an extreme Kerr black hole rjj — Tg. 
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The static limit ve is situated at 

rs = - fl + Vl-a^ cos2 d\ (5) 

where 9 is the angle to the polar axis. The region between rn < r < te is called 
the ergosphere. Within the ergosphere (which is ellipsoidal in shape), space-time is 
dragged along in the direction of the hole's rotation (frame dragging), so that no 
static observer can exist and a particle has to co-rotate with the hole. The angular 
velocity J7jy of a black hole is defined as the angular velocity of the dragging of 
inertial frames at the horizon and given as 

0„^«(^). (6) 

SO that for an extreme Kerr black hole fin — c/2rg. 

A rotating black hole possesses an additional source of energy, which for low spin 
parameters a is given by Srot — (l/2)/_f/^|f = {1/8)0^1140^ where the black hole's 
moment of inertia Ih = Mrjj has been employed. For general spins, the maximum 
rotational energy i^iot of a rotating black hole can be expressed as the difference 
between the total mass-energy Mc^ and its irreducible mass Mi„c^ , see refsEHMl, 

E,ot = {M - Mi,,)c^ = Mc^ (l - yWo-M/r^^o^) < 0.29A/c2 . (7) 

Thus, if supermassive black holes would indeed be rapidly rotating and if they 
would be able to dissipate all of their rotational energy over a characteristic life 



time iEdd = 4.5 x 10^ yr (see Eq. 11 ), this could yield a considerable luminosity 
output 

L„. = |=i<3.6xlO«(jjiL-)e,g/., (8) 

comparable in magnitude to the Eddington luminosity (see Eq. ^). The maximum 
output is, however, very sensitive to the spin parameter a. For moderate spins, e.g., 
a — 0.5, it is already about an order of magnitude smaller (-Eiot = 0.034Mc^), see 
Fig. ^ for illustration. 

Note that the most general, stationary black hole metric is the Kerr-Newman metric, 
where the black hole is described by its mass M, electric charge Q and spin angular 
momentum J. Special cases are the Kerr (Q = 0, J 7^ 0), the Reissner-Nordstrom 
( J = 0, Q ^ 0) and the Schwarzschild {J = Q = 0) black hole. Charged black holes 
are commonly thought not to be of astrophysical importance (but see also, e.g., 
refsEIE3), as the surrounding plasma is usually expected to rapidly neutralize any 
charge imbalance (e.g., refsESElljr] 



^The metric function implies the constraint {J / M)"^ -\-{GQ'^ / c^) < (GM/c)'^, so that the maximum, 
theoretically possible charge would be Qm = G^'^M. However, real astrophysical black holes are 
expected to possess only negligible charge (i.e., Q ~ GMmp/e ^ 10~^^Qm) as otherwise electric 
forces would dominate over gravitational ones, enabling charge separation and attracting opposite 
charge to neutralize the imbalance. This implies that the magnetic field produced by a charged 
black hole is Umited to B ^ mpc'^/{erg) ~ 2 X lO-'^(lO**Af0/M) G. 
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Fig. 3. Maximum cxtractablc energy of rotating black liole as a function of the dimcnsionless 
spin parameter a. 



2.2.2. On the Distribution of Black Hole Spins 

Despite their obvious significance, the spins of astrophysical (supermassive) black 
holes have been proven difficult to pin down unambiguously. This may appear some- 
what unfortunate, as it e.g. directly affects inferences about the maximum rota- 
tional energy that could possibly be extracted by some electrodynamic processes 
(cf. Eq. [?]). Naively, one expects that supermassive black holes can be spun- up 
by prolonged accretion of angular momentum or by the merger of black holes with 
similar masses (having high orbital angular momentum). On the other hand, if a 
rotating black hole would be able to emit a strong (persistent) spin-powered jet 
of, e.g., Lj ~ lO'*^(Af/lO*M0) erg/s, the hole's rotation could spin down on an 
e-folding timescale E^ot/Lj ^ 10^ yr comparable to the typical AGN lifetime. Then 
again, prolonged accretion or further merger events could help to sustain high black 
hole spins. Prolonged accretion is expected to lead to high spins (a ~ 1), 1251 whereas 
major mergers tend to produce average spins of (a) ~ O.?!'^^' Provided accretion is 
not inherently random and episodic, cosmological N-body simulations suggest that 
the final black hole spin depends almost exclusively on the accretion history and 
very little on the merger history.!^ Accordingly, more massive black holes (i.e., those 
hosted by elliptical galaxies) should be characterized by higher spin parameters (av- 
erage (a) > 0.8 at z = 0) compared to those residing in, e.g. spirals. However, the 
underlying assumption of a spin evolution via prolonged accretion has been ques- 
tioned recently.^® In the standard picture, Lense-Thirring precession is assumed to 
always ensure co-alignment of the disk with the black hole spin, and thus to lead 
to spin-up via prograde accretion.'22] However, this needs not necessarily to be the 
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case if the disk is sufficiently small. Instead, accretion may then be random and 
episodic in character, so that the spins are expected to readily adjust to average 
values (a) ~ 0.1 — 0.3. Such a spin-down scenario seems in fact supported by recent 
studies of the radiative efficiency from accretion onto supermassive black holesrl 
(i) Existing quasar survey data suggest a significant decrease of the radiative ef- 
ficiency from high {rj ^ 0.2 — 0.3 at z ~ 2) to low redshift (ry ~ 0.03 at z = 0). 
While the spin evolution at z ^ 2 may be driven by major merger events, the de- 
crease of the radiative efhciency towards low z seems to support an interpretation 
according to which black holes are fed by accretion events that tend to spin down 
the black hole over time (a ^ at z = 0).'22 (ii) Similarly, a rather low mean efh- 
ciency {rj) ~ 0.07 has been inferred based on the cumulative energy density emitted 
by AGNs over the age of the Universe. The results indicate (a) ~ (0.2 — 0.6) and 
suggest that supermassive black holes are on average not rapidly spinning during 
accretion.133 Note that if Maxwell stresses do (contrary to what is usually assumed) 
not vanish at r-^g, as in fact seen in recent MHD disk simulations,'^ this could 
lead to an excess dissipation, lowering the above inferred spins even more (possible 
differences seem to be small, though, at least in the case of a thin dislP^. 
Yet, while on average rather moderate values of a seem to be preferred, a sub- 
population of AGNs with high spins may well exist. In fact. X-ray observations 
of relativistically broad iron Ka lines from the black hole's vicinity, e.g., in the 
prominent Seyfert 1 galaxy MCG-6-30-15, for which a > 0.9 has been inferred, 
indicate that at least some radio-quiet (!) AGN may harbor rapidly spinning black 
holes (for review, see refsji^'^ '^^)F]lf indeed the case, then this implies that a high 
black hole spin is not a sufficient condition for radio-loudness as initially suggested 
by Wilson & ColberlP2l. On the other hand, evidence for rather low black hole 
spins in the galactic microquasar sources XTE J1550-564 (a ~ 0.1) and A0620-00 
(a ~ 0.3) also challenge the common perception that black hole spins and relativistic 
jets are necessarily connected.HSl Fender et al.H^have recently concluded that at least 
for galactic black hole X-ray binaries (which are often believed to be scaled-down 
versions of AGNs), there is no evidence for either the jet power or the jet speed 
being related to the spin of the black hole. If this is indeed the case, it may point 
to the fact that the jets rather emerge as centrifugally-driven (Blandford-Payne) 
outflows from the accretion disk. 
With this caveat in mind, estimates for the jet kinetic power (e.g., from intraclus- 



''The radiative efficiency r] is usually defined as 77 = 1 — ems/c^, where ems = c^(l — rs/Srms)^''^ 
is the binding or specific energy of a test particle rotating in the innermost stable circular orbit 
at rms, and becomes rj = 0.057 for a Schwarzschild and rj = 0.42 for an extreme Kerr black hole 
(a = 1). Note that while accretion by itself would spin up the black hole rather quickly to its 
maximum value a = 1, photon capture will prevent complete spin-up and limit the maximum spin 
to a = 0.998. While this is still very close to a = 1, it has a significant effect on the maximum 
possible radiative efficiency, reducing it to t; = 0.3.1221 

•^Note however, that recent modelling of the broad iron lines in the Seyfert- 1 AGNs Fairall 9 and 
Swift J2127. 4-1-5654 seems t o favo r more moderate spin values a ~ 0.6, making MCG-6-30-15 
appear somewhat exceptionalEZEH 
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ter X-ray cavities inflated by the radio jets) and for the central black hole mass 
have been employed to constrain the spin parameters for selected samples of radio 
galaxies. Some of these studies conclude that the radio-loud dichotomy of FR I and 
FR II radio galaxies can be reproduced by a transition in accretion mode from a 
standard disk to an ADAF, with the black holes for both classes rapidly rotating 
(a ~ 0.9) HEU Others instead find that the FR I type and FR II sources they studied 
are consistent with spin values of a ~ (0.01 — 0.4) and a ~ (0.2 — 1), respectively, 
and report evidence for a decrease of the black hole spin in FR II with decreasing 
redshift."^ At present, further studies are certainly needed to clarify this issue. 

2.3. Accretion onto Black Holes 

2.3.1. Spherical accretion and Eddington constraints 

The luminosity of an accretion-powered object of mass M cannot grow indefinitely. 
In the case of steady spherical accretion, the maximum luminosity is limited to the 
Eddington value LEdd given by the balance between (outward-directed) radiation- 
pressure and (inward-directed) gravitational force, 

^Bdd = -^^ = 1.25 X 10- [^-^ erg/s . (9) 

The critical accretion rate required to sustain the Eddington luminosity, assuming 
a typical 77 = 0.1 efficiency of conversion of mass to radiant energy, is then defined 
as 



The Eddington timescale (sometimes called the Salpeter timescale) is the e-folding 
growth time for the mass of a black hole accreting at the Eddington rate 

tEdd = -^-4.5xl0V- (11) 

The fiducial Eddington magnetic field upper limit i?Edd h^s an energy density 
equal to the radiative energy density at rg of a body emitting at ^Eddj i-e., 
B = (2LEddAg2c)V2, 



BEdd ^ 6.1 X 10" ( ^^ ) G. (12) 



2.3.2. Standard, steady state accretion flows 

In general, the structure of an accretion flow is determined by the balance between 
gravitational heating and cooling. Its concrete structure therefore depends on what 
kind of heating and cooling processes are assumed to be dominating. In the classical 
approach of Shakura & Sunyaev - , commonly referred to as the standard accretion 
disk model, the disk is flat (geometrically thin: H -^ r; H the half thickness of the 
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disk, r the radial distance) and opaque (optically thick in vertical direction). Viscous 
(frictional) stresses^ are assumed to convert gravitational potential energy into heat 
that is released locally in the form of a thermal black body spectrum B^{Tcs[r]) 
(i.e., possible modifications by advection or a jet are neglected) with energy flux per 
unit surface area of F{r) = tt J B^{Tcs[r])di' = <TTcs{r)'^ and characteristic effective 
temperature 



Teff(r) 




6.3 X 10 



Far from the inner edge, r ^ r^, the temperature approximately obeys reff('') oc 
j,-3/4 ^g gravity is stronger on smaller scales, the disk surface is hotter in the inner 
region. The radially integrated disk luminosity is 



id 



/' 

•'ri. 



2-KrF[r) dr 



GMM 
2nn 



Mc^ 



(14) 



Half of the total gravitational energy release GMM/rin is thus radiated away, the 
remaining half being accounted for by the kinetic energy of the orbiting material. 
The emerging disk spectrum F^, ex J Bi^(Tcff[r])rdr can be obtained by integrating 
Bij over the surface of the disk. The resultant spectrum is shown in Fig. HI It initially 
rises with i/^, then flattens (provided the disk is large enough) in the intermediate 
region to v^^^, before it decays exponentially towards the highest frequencies. For 
characteristic AGN parameters, the emission from the inner part is maximized at a 
frequency 



2.82kTcS 



5.9 X 10 



15 



105 K 



Hz. 



(15) 



suggesting that this type of accretion flow could lead to a 'big blue bump' feature, 
i.e. a signiflcant optical/UV disk contribution to the observed SED. 
As a consequence of the viscous interactions between adjacent layers in the disk, 
angular momentum is transferred outwards, while mass falls (accretes) inwards. 

2.3.3. Inner edge of the accretion disk 

In the classical picture, the radius of the inner edge r-^ of the standard disk extends 
down to the (innermost) marginally stable circular orbit r\a — ?'ms- Inside r.^^ 
space-time gets so strongly curved that stable orbits no longer exist. A particle 
thus unavoidably falls into the black hole. This is why the region inside r^s is often 



"With the r(/)-coinponcnt of the shear stress tensor taken to be proportional to the total (gas and 
radiation) pressure p, i.e., t^^ = ap, a < 1 (a-prescription). 
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hv/kT(ro„,) 

Fig. 4. The characteristic continuum spectrum of a standard, geometrically thin, optically thick 
disk for two different outer radii rout- The spectrum is normalized to hu/kT{rout)- If the disk is 
not large enough, e.g., for rout/^in = 100, the i^^'^— power law part may not become apparent. 

called the "plunging region" as gas can spiral into the black hole without further 
loss of angular momentum. The analytical solution for the innermost stable orbit 
as a function of the dimensionless spin parameter a is given bjSl 



3 + Z2T a/(3-^i)(3 + Zi+2Z2) 



(16) 



where 



Zi = l + il-ay^^ (l + ay + il-a) 



a/3 



Nl/3 



3a2 



zl 



Hence, for a Schwarzschild black hole (a — 0), 



extreme Kerr case (a = 1), 



I = 3rs = 6rg 
rs/2 — Tg (prograde orbit) or 



(17) 
while for the 



9rg (retrograde orbit). As mentioned above, the disk's radiative efficiency is strongly 
dependent on the position of its innermost radius. 



2.3.4. Radiatively inefficient accretion flows 

In the standard approach, cooling is so efficient that all the energy released through 
viscosity is radiated away locally. This may not always be the case and in fact 
we know steady disk solutions where this assumption is violated (for a review, see 
Kato et al.'^Zl). In the case of a radiatively inefficient accretion flow (RIAF), for 
example, most of the energy released is stored within the flow and transported 
inward with accretion. While the disk thus locally undergoes advective cooling, the 
flow itself becomes very hot. When the density of the accretion flow falls below 
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a critical value, the flow becomes optically thin and the efficiency for radiative 
cooling (e.g., bremsstrahlung, which dominates at non-relativistic temperatures) 
becomes small. The classical RIAF prototype is the two-temperature, optically-thin 
advection-dominated accretion flow (ADAF) model (for review, see refsEHMl). By 
assumption, the released viscous energy is assumed to go primarily into ion heating, 
while cooling is mainly by the electrons. If the densities are sufficiently low. Coulomb 
coupling between protons and electron is weak and the viscous energy transfer from 
protons to electrons becomes very small, limiting the amount of energy that can be 
lost by radiation. This introduces a critical accretion rate 

M ~ a^A/Edd (18) 

above which an ADAF cannot exist (typically, a ^ 0.2 — 0.3). Because of the very 
high temperature {Tp ~ GMmp/Sksr - 3 x 10^'^ (rg/r) K for ions, Tg - 5 x 10^ K 
for electrons) ADAFs are marginally geometrically thick with scale height H/r ^ 1 
at every radius, independent of the accretion rate. The total radiative luminosity of 
an ADAF is found to be proportional to the square of the mass accretion rate, in 



contrast to the standard thin disk where L ex M (eq. 14 ), i.e. 

Ladaf - ^^^ (^^] Mc' oc M^ (19) 

(Mahadevan 1997), and thus can be much smaller than the luminosity expected 
from a standard disk. The broadband spectrum emerging from an ADAF can cover 
the energy range from radio to hard X-ray frequencies, see Fig. [51 The radio to sub- 
mm regime is typically produced by synchrotron emission of relativistic thermal 
electrons and in its optically thin part rises with L^ ex v'^^^ W^ The emission at the 
highest peak frequency Vp comes from the innermost part of the disUSHEEl 



..W.IO'Mt^J UIm^ y^} (-) H., ,20) 



r 



while the lower frequency synchrotron emission originates from further out. Comp- 
ton upscattering of the low-energy synchrotron photons by the relativistic ther- 
mal disk electrons results in a hard power law tail extending up to energies of 
hv ~ fcTe ~ (100 — 500) keV. Compton scattering becomes less important with de- 
creasing accretion rates, and the spectrum becomes steeper (softer). For sufficiently 
low M distinct Compton peaks can appear. This is related to the increase in elec- 
tron temperature with decreasing accretion rates. For low M, the X-ray spectrum 
is dominated by bremsstrahlung emission (with a non-negligible contribution from 
large radii) due to electron-electron and electron-ion interactions, again up to a 
maximum energy of ~ fcTg, beyond which the spectrum falls off exponentially.ED 
A number of additional scenarios have been discussed in the literature that extend 
and complement the classical ADAF picture described above. One important mod- 
ification concerns, e.g., jets and winds. Since for self-similar ADAFs the net energy 
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Fig. 5. Sketch of a typical ADAF spectrum around a supermassive black hole. Synchrotron emis- 
sion (S) of the relativistic thermal electrons produces a peak at around lO^'^ Hz. At moderate 
accretion rates, Comptonization results in a power-law tail (C) extending into the X-ray domain. 
For low accretion rates, distinct Compton (Cl=once Compton scattered, C2=twice Compton scat- 
tered) peaks appear. The X-ray/soft gamma-ray regime is usually dominated by bremsstrahlung 
emission (B), peaking at hu ~ fcTe, and with luminosity scaling as L^r oc m^. 

of the accretion flow can remain positive and unbound to the central object (posi- 
tive Bernoulh number),'^' strong outflows may well form and modify the emergent 
spectra, cf. the ADIOS model by Blandford & BegelmarPSl. While the positiveness 
of the Bernoulli number is a necessary (but not sufHcient) condition for outflows, 
2-D simulations indeed suggest that for large viscosity parameters (a ~ 0.3) strong 
outflows can occur.^li5.^ If mass is supplied by an outer standard disk, ADlOS-like 
outflows may however only occur for transition radii rt ~ lOOfs- 



2.4. Magnetic Field Strengths 

In the standard disk case, we can derive a characteristic (maximum local equipar- 
tition) magnetic field strength close to 
that of the radiation (i.e., B^ 



Tg, whose energy density is comparable to 



in = aT^g/c) which gives (cf. eq. 13 ) 



Br 



4.4 X 10^ 




(21) 



For M — AfEdd this is of about the same strength as the Eddington equipartition 
magnetic field -BEddi cq. (12 1, obtained by assuming that all of Lgdd emerges from 



In the case of an advection-dominated accretion flow (ADAF) on the other 
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hand, the characteristic equipartition magnetic field is given by B'^/Stt ~ 0.5pc^. 
Here Cg denotes the isothermal sound speed, p oc M /{r'^Vr) is the density of accreted 
matter, Vr ~ Q.havf is the typical radial infall speed and Vf — (GM/rY^^ the free 
fall velocity. When scaled to the AGN mass range, this become^^ 



As an ADAF cannot exist above AI ^ a^AfEdd, cf. cq. (18 1, the maximum Ba is 
comparable to the maximum Br in the standard case. 

Note that the poloidal field strength in the magnetosphere may be significantly 
smaller than these disk field values. Taking Br to be characteristic for the poloidal 
magnetospheric field component, for example, may then well overestimate the pos- 
sible power of a Blandford-Znajek-type process. It has in fact been argued,!^ that 
the large-scale poloidal field Bpd threading the disk should be related to the small 
scale field produced by dynamo processes via Bpd ^ (-ff/'')-Sdynamo, where H/r is 
the typical disk scale height, which for a standard disk is much smaller than unity 
(cf. however also ref.'^. In any case, as the field threading the hole is generated by 
currents in the disk, it seems unlikely that its strength should significantly exceed 
the field threading the inner disk. 

2.5. Black Hole-Jet-Magnetospheres 

Once a spinning black hole or disk becomes threatened by an ordered (large-scale 
poloidal) magnetic field, it can build up a rotating magnetosphere filled up with 
currents. The required magnetic fiux could be advected from the interstellar mag- 
netic field via the accretion process and/or be produced (amplified) through dynamo 
actions in the inner accretion diskF] While this picture may appear qualitatively ev- 
ident, the detailed electromagnetic structure of the magnetosphere is a highly com- 
plex problem. Most research has been infiuenced by the seminal papers of Blandford 
& Znajek (BZ)P21 and Blandford & Payne^. Roughly speaking, the main difference 
concerns the question whether the source of rotational energy and angular momen- 
tum lies with the black hole's rotation or with the accretion disk. Methodologically, 
this often translates into two different approaches (for review and further references, 
see e.g. refsHUS): 

In the first, the force- free electrodynamic (FFE) approximation, the magnetospheric 
structure is assumed to be governed by Ampere's and Faraday's law and to be 
force-free everywhere (i.e., dominated only by magnetic and electric stresses). Ac- 
cordingly, gravitational and inertial forces are neglected (massless limit of MHD) 



•^An alternative scenario relates to the Poynting-Robertson radiation force acting predominantly 
on the electrons in the accretion disk. As a consequence, electrons lag behind the protons, and a 
toroidal electric current develops that can become sufficiently large to produce significant poloidal 
magnetic fields.ESl 
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and the current density j perpendicular to the local magnetic field is determined 
from the force-free condition pEj_ + j x B/c = 0, where p is the charge density!^ 
This approximation quite simplifies the standard problem in that it allows to solve 
for the electromagnetic field structure without having to solve for the plasma dy- 
namics. Magnetic field lines may then be thought as (quasi-rigidly) rotating with 
the angular velocity J7 of the horizon or the inner disk (Ferraro's law). Approaching 
the light surface where |f7 x t^ = c the field is swept backwards opposite to the 
sense of rotation and a toroidal twist is introduced. By assumption, plasma in the 
magnetosphere is streaming outward at almost the speed of light, its only role being 
to provide the currents required to sustain the magnetic field. 
In the second approximation, the relativistic MHD (single fluid) approach, the force- 
free condition is modified to allow for inertial and pressure terms in the equation of 
motion. If plasma loading is not negligible, as may be the case for field lines connect- 
ing to the disk, such a modification appears unavoidable. The system is then closed 
by Ohm's law asserting that the current density j is proportional to the electric field 
E' = E+vx B/c seen in the frame where the fiuid is at rest. In the high conductivity 
(ideal MHD) limit, Ce — > oo, Ohm's law implies j/cfe = E + vx B/c — > where wis 
the fluid velocity. Within the Alfven surfacq^ plasma that becomes attached to the 
field lines behaves like a bead on a wire and is fiung out by the centrifugal force given 
suitable field line orientations (e.g., inclinations wrt to disk launching of less than 
60° in Newtonian approximation, and nearly 90° in the extreme Kerr limit EMMIl). 
For stationary and axisymmetric configurations, the equation of motion (which in 
the low inertia limit just reduces to the force-free condition) projected perpendicular 
to the field lines gives the famous Grad-Shafranov equilibrium (ideal MHD) equa- 
tion. The so-derived Grad-Shafranov equation (GSE) is a highly non-linear partial 
differential equation and (while simply re-expressing the cross-field force balance) 
highlights that the structure of the magnetosphere essentially depends on the as- 
sumed field line rotation and the current distribution in the magnetosphere EMEl 
In the low inertia (high magnetization) limit and flat Minkowski space, the GSE 
(without differential rotation) reduces to the original pulsar equation and describes 
the solution of the time- independent set of FFE equations. An example of a possible 
magnetic field structure emerging from the disk is shown in Fig. [61 
Both, the FFE and the MHD fluid approach, essentially rely on the assumption that 
charge separation is/remains negligible everywhere such that E ■ B ~ Q is satisfied 
over the length scale of interest. This could be a reasonable assumption for e.g., 
field lines connecting to the disk and the surrounding plasma: 
In the laboratory frame, the rotation of the magnetic field will induce an electric 
field E = —Viot/c X B. If not screened, this could in principle be tapped for the 



Defined in the non-relativistic case as the surface where the plasma kinetic energy density equals 
the magnetic field energy density. Relativistic flows remain Poynting-dominated at the Alfven 
surface, with the Alfven surface only being related to the propagation of Alfven waves, that 
become static at this location. 
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Fig. 6. Poloidal magnetic field structure of a coUimating, relativistic jet, obtained from the numer- 
ical solution of the stationary-state, axisymmetric relativistic MHD equations. The magnetosphere 
is assumed to be force-free. Dipolar magnetic flux carried by the accretion disk is stretched by 
a disk wind into the upper hemisphere. Left: Assuming weak differential (quasi-rigid) rotation of 
foot points of the field lines. Right: For strong diff'erential rotation. Rq denotes the asymptotic 
light cylinder radius. From Fendt and Memola^. 

acceleration of particles. By Gauss' law, the induced electric field is supported by 
a local charge density pe — V ■ E /At:, corresponding to a particle number density 
(commonly referred to as the Goldreich- Julian [GJ] density^^) well inside tl of 
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where fi = c/rL, tl is the typical light cylinder scale with tl = rjTs, and 7y ~ (1 — 10) 
for disk and -q — 1/a ^ 1 for (extreme) Kerr black hole magnetospheres. Here, 9 
denotes the field inclination. If the environment is plasma-rich {n > uqj), the ability 
of charges to move freely along the field lines will ensure (at least inside tl) that 
E-B — >■ 0, i.e., the parallel electric field component essentially becomes screened. (As 
a natural consequence, gap-type particle acceleration in the induced electric field 
will be suppressed.) Applicability of the MHD approximation thus basically requires 
that the charge density is always larger than the critical Goldreich- Julian density. 
For self-consistency one may thus always wish to verify that this is satisfied within 
a numerically calculated flow.'^ The inner regions of AGNs are usually sufficiently 
plasma-rich to make this a useful working assumption for disk-driven outflows. If 
one conservatively assumes, for example, that accretion occurs close to the free 
fall timescale, i.e., t^ ^ r/vr, Vj. ~ Q^^/j ^/ — (GM/r)^''^, with a the viscosity 
parameter, the characteristic particle number density n ^ Mta/TTr'^nip close to r^ 
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iSS disk 




Fig. 7. Schematic drawing of a possible field structure around a rotating black hole, with a black 
hole-driven wind surrounded by an outflow from the inner parts of the disk. The accretion flow 
may switch its configuration to a standard disk at some transition radius rj. 
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i.e., orders of magnitude higher than the minimum GJ density required to screen 
the fields. The physical situation is less obvious for black hole-driven jets. Here, 
the nature of the plasma source that replenishes charges which escape the system 
along open magnetic field lines (both, outgoing to infinity and ingoing across the 
horizon) is less well understood, and charge-starved funnel regions or gaps can occur 
for which the ideal MHD description may not be adequateE3Zl Plasma injection is 
then most likely be related to electron-positron pair production (cross-field diffu- 
sion being negligible), either due to direct 77-interactions of MeV photons in a hot 
accretion flow and/or initiated by, e.g., the Compton up-scattering of ambient pho- 
tons by electrons/positrons accelerating in a vacuum gap (e.g.,refs P° l ^^ l ^^ l ^^ l ^^ l ^^ l ^ l). 
For a number of conditions, this is believed to ensure a plentiful supply of charges, 



cf. also eq. (68), so that the plasma around (e.g., above the gap) a "typical" AGN 
is commonly also considered to satisfy quasi-neutrality (n'^ + n~ 3> n" — n'^)?^ 



Global 3d-relativistic MHD simulations of black hole accretion disks and mag- 
netospheres have in recent years benefited from a proliferation of numerical ap- 
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proachesp] The results are widely considered to answer, at least partially, some 
of the longstanding questions with respect to magnetically-dominated (Poynting- 
flux-dominatedpiJ and matter-dominated jets:-''' Accordingly, the simulations con- 
firm the general picture of a two-component sheath jet configuration with a fast, 
magnetically-dominated inner jet being surrounded by a slower moving, matter- 
dominated jet. As it appears, relativistic magnetically-dominated jets may in fact 
need an external agency (e.g. a disk wind) to facilitate efficient coUimation (cf. also 
refsE3121). There are indications that the jet formation process does not require or- 
dered, large-scale magnetic fields to be fed in from larger disc scales,*^ and that the 
resultant outward-directed electromagnetic energy flux is consistent with the flux 
being powered by the Blandford-Znajek process. ^^ The first point, however, may 
need some further confirmation as other simulations still seem to need large-scale 
seed fields (e.g., ref.^^). 

For the obtained matter-dominated jets, on the other hand, the main driving force 
seems to be a high-pressure corona (no evidence for the operation of a Blandford 
& Payne-type process has been found yet), where strong (gas plus magnetic) pres- 
sure gradients ensure that matter is pushed into the jet (leading to extended mass 
loading) and accelerated along it.l^This comes along with relatively slow jet speeds 
{Tj ~ 1.5) and large opening angles. 

While very instructive, these global simulation results still face (apart from the usual 
methodological constraints, see e.g., ref.l^ a number of limitations whose infiuences 
need to be carefully assessed in order to allow detailed contact with observations. 
Among others, this relates to the facts that: (i) Not all of the global numerical 
schemes employed are fully conservative (e.g., in Hawley et al. only the internal, 
but not the total energy is explicitly conserved), (ii) Radiation fields, although ob- 
servationally important, are as yet not fully included (e.g., in McKinney). This 
could obviously modify the fiow dynamics and considerably influence the heating 
and cooling of the plasma. In fact, it has been argued that for powerful quasar- type 
sources, radiative driving may play an important role in shaping the outflow.l^Zl 
(iii) Usually, the plasma content is assumed to be conserved, which may not be 
appropriate in the vicinity of the black hole, where pair creation is occurring all 
the time (e.g., refsEMl). (iv) The results can depend strongly on the chosen ini- 
tial and boundary conditions. Since magnetic fields are, e.g., divergence- free, they 
seem never able to fully forget their (artificial) initial configurations, making the 
resultant jets very sensitive to these (observationally poorly known) conditions.'^ 
(v) While global approaches rely on the ideal (infinite conductivity/zero resistivity) 
MHD limit, longterm simulation results seem to provide evidence for (artificial) 
reconnection-driven magnetic island formation. The occurrence of magnetic recon- 
nection could significantly change the global field configuration and the plasma 
dynamics®^ (for an illustration of the role of finite resistivity in the Newtonian 



^The interested reader is referred to the relevant reviews by, e.g. Fragile^ and Krolik & Hawley^ 
for more details and references. 
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limit, see refpfl for the special relativistic case see refpU). It is not clear to which 
extent this might invalidate some of the numerical results obtained so far with ideal 
MHD. 

Nevertheless, for the purpose of our investigation here, we may visualize a sim- 
plified field structure as shown in Fig. [7| 

3. PARTICLE ACCELERATION IN BLACK HOLE-JET 
MAGNETOSPHERES 

3.1. Direct electric field acceleration and the Blandford- Znajek 
process 

A black hole embedded in an external poloidal field Bp and rotating with angu- 
lar velocity fin — a(c/2r//), Eq. (p|, will induce an electric field of magnitude 
\E\ ^ {^HfH)Bp/c corresponding to a permanent voltage drop across the horizon 
of magnitude <& ^ rH\E\ = {a/2)rHBp. In terms of the electric circuit analogy, a 
black hole can act like a unipolar inductor (battery) with non-zero resistance, so 
that power can be extracted by currents flowing between its equator and poles. For 
reasonable astrophysical parameters, the potential drop (noting that one statvolt 
[cgs] corresponds to 300 V [SI] ) is of order (cf . ref.I^ 

*~2xlO»„(l + vT^)(jJ^)(j^)|Vl (2.) 

If a charged particle can fully tap this potential, acceleration up to ultra-high en- 
ergies oi E = Ze^ - 3 X m^^Z{M/10^MQ){Bp/lQ'^ G) eV may become possible. 
Rapidly spinning and massive black holes have thus often been considered as as- 



trophysical candidates for the energization of UHE cosmic rays (cf. § 5.2.2). As 
noted above, however, the charge density in the vicinity of accreting black holes 
may well be so high that a significant fraction of this potential is screened and thus 
no longer available for particle acceleration (see, however, also ref.l^. It seems thus 
more appropriate to define an effective potential where the available length scale 
(gap height h) is explicitly taken into account, i.e.. 



*e-^//(^-j ri,Bplc=^-j $. (26) 

Accordingly, the characteristic rate of energy gain can be expressed as 

dE 

— = Ze<^ec/h , (27) 

at 

with corresponding acceleration timescale tacc — jrnQch/lZe^e)- 
Equation ( [25| yields a strong electromotive force (EMF) which in the Blandford- 
Znajek (BZ) scenario is assumed to drive electric currents along the fields and to 
generate a power output Pbz = ^I — ^'^ /{^Zs). As the event horizon docs not have 
perfect conductivity but instead acts like free space with respect to the propagation 
of electromagnetic waves, its surface resistivity Zs corresponds to the impedance of 
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free space (377 Ohm — Att/c [stat-ohm], with c in cgs)lH2E2lThus, if current closure 
can be achieved in the black hole environment (for field lines emerging between 
pole and equator, i.e. with resistance AZg ^ Zg/A ^ 100 Ohm), a considerable 
electromagnetic (Poynting-flux-driven) power output of the order 

PBZ^5xlO«a2(l + v/rr^ffc('-il-V('-§^V erg/s. (28) 



might be generated, where a is the dimcnsionless spin parameter and A; ~ 1 is 
a numerical factor accounting for the field geometry. To support such an elec- 
tromagnetic power output, however, an electric current / = (Pbz/^-^s)^ ^ 
7 X lO^^statampere ~ 2 x lO^^A must fiow within the magnetosphere, requiring 
a minimum steady pair injection (creation) rate of dNf,/dt = I/Q ~ 10^^ par- 
ticles/s. The implied (minimum) particle number density would then be close to 



the Goldreich- Julian density estimated in eq. ( 23 1 . Even if one neglects cross-field 



diffusion, this density could, as shown later on, most likely be maintained by pair 



creation alone, cf. eq. (68) 



The above analysis, when put on a more rigorous basis (cf. the "Membrane 
paradigm'^f^Sl) ^ reveals a number of subtleties that have been the center of much 
attention and debate. To mention a few: 

(1) In reality, the BZ power depends on the assumed magnetic field topology. The 
angular velocity of the field lines Qf, for example, does not simply coincide with 
the angular velocity of the black hole Hh, but is instead dependent on the global 
field structure. Even if the impedance-matched case (where the resistance of the BH 
"battery" equals the one of the "load" at large distances, i.e., where flf — 51 ///2) 
becomes established with time such that maximum power extraction is ensured,!^ 
the numerical factor k might well be smaller than one (e.g., k < 0.2 for a homoge- 
neous magnetic field, see ref.l^. In general, for an Eddington-limited system, the 
BZ process cannot extract energy at a rate much larger than 0.29 iEdd/^, where 
77 ~ 0.1 in the canonical conversion efficiency (cf. eq. (sl; also refslHMZI). Recent 
studies of powerful FR II radio galaxies suggest that these sources can have jet 
powers comparable to the accretion power P ~ AIc^ (e.g., ref.^-). This is much 
higher than expected from current, global relativistic MHD simulations of the BZ- 
process, that start off from a no-net-magnetic flux conflguration, and may indicate 
that large-scale magnetic flux needs to be accreted in order for the process to oper- 
ate efficiently enough (see rei^ioi details). 

(2) The "membrane" analogy, according to which the black hole horizon can be 
regarded as unipolar inductor with non-zero resistivity, has been fundamentally 
criticized on physical grounda^^ni (cf. also ref .li^U for a response). As the horizon is 
causally disconnected from the external space and the outgoing flow, it actually 
cannot act as a unipolar inductor, which has led to general doubts concerning the 
nature of the BZ process. There are at least two (related) responses to this prob- 
lem: (i) First, within the Grad-Shafranov MHD approach it can be shown that the 
" boundary condition" determining the potential energy loss of a rotating black hole 
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is in fact dictated by the physical parameters in the pair creation region outside 
and not at the horizon IMli^^ This gives about the same power output and iUustrates 
that the nature of the BZ process does not have to rely on a questionable causal 
connection between the event horizon and the outer magnetosphere. (ii) Secondly, 
the physical source of the BZ process does not - when correctly interpreted - lie 
with the horizon itself, but with the ergosphere. As plasma is unavoidable forced 
into co-rotation within the ergosphere, the magnetic field (assumed to be frozen in 
the plasma) is forced into co-rotation as well. While the horizon is indeed causally 
disconnected from the outgoing wind, the ergosphere is not and the BZ mechanism 
can still be recovered ISIEMI 

(3) In the BZ scenario, black holes are assumed to be surrounded by plenty of 
charges, ensuring that an electric current can flow along the poloidal magnetic 
fields. This poloidal current then results from the screening of the vacuum electric 
field. On the other hand, to sustain a stationary current, the electric field should 
retain small unscreened components and not become completely screened by a re- 
distribution of electric charges. Otherwise, the magnetospheres would no longer be 
able to drive electric currents (cf. ref.^'^°). As it turns out, however, such a final state 
appears not possible within the ergosphere of a rotating black holeP^H The electric 
field, which is induced "gravitationally", cannot become completely screened by 
any static distribution of electric charges, i.e., rotating black holes do not allow 
stationary solutions with completely screened electric fields and vanishing poloidal 
currents. 

(4) The BZ process of extracting the rotational energy of black holes appears to 
be very similar, but not simply equivalent to the mechanical Penrose process.E^Sl 
Because plasma entering the ergosphere is forced to co-rotate with the black hole, 
magnetic field lines penetrating the ergosphere will get twisted by the differential ro- 
tation of the plasma. This twist of magnetic field lines propagates outwards along the 
field lines (against the infalling plasma flow) as a torsional Alfven wave, establishing 
an outgoing electromagnetic energy flux. The magnetic tension of the bent field lines, 
on the other hand, decelerates the plasma in the ergosphere, imparting opposite an- 
gular momentum. The ergospheric plasma thus moves to lower energy orbits and 
the total (hydro- and electrodynamic) energy-at-infinity quickly decreases to nega- 
tive values.'^^ While it seems that initially the resultant Penrose energy extraction 
process essentially operates via negative hydrodynamic energy-at-infinity, longterm 
simulations suggest that regions of negative hydrodynamic energy-at-infinity even- 
tually disappear, with the electromagnetic BZ process still continuing.lJ^ A more 
general analysis suggests, however, that the condition, for which the direction of the 
Poynting energy fiux (at infinity) is opposite to that of the (local) field flow, is not 
simply equivalent to the Penrose one for negative energy-at-infinity (see, e.g. ^^^ for 
more details). 

(5) According to recent results based on general relativistic MHD simulations, the 
maximum BZ jet power appears to be sensitive to the assumed jet-disk geometry. 
In the presence of, e.g., a thick disk (with effective thickness (H/R) ^ 1) the jet 
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power may actually scale more strongly with black hole spin, i.e., Pbz oc a^ instead 
of Pbz oc a^ for a thin disk.^''^ If this is confirmed by further research, it may offer 
important clues for understanding the radio-loud/radio-quiet dichotomy in AGN. 

3.2. Direct electric field acceleration close to the light cylinder 

In early applications of the BZ process, charged particles (electric currents) are 
considered being able to cross the field lines only in the "load" region far away from 
the black holcESl Accordingly one may expect significant particle acceleration to 
occur mainly within that region. However, this needs not necessarily to be the case 
as has been demonstrated for pulsar magnetospherespii' Within ideal MHD, the 
structure of the magnetosphere is essentially determined by the supposed field line 



rotation and the assumed electric current distribution (cf. § 2.5). If the longitudinal 
electric current would become sufficiently small, the freezing-in condition will break 
down so that electrons, deviating from the field line, may become able to sample 
the full induced electric field and thereby get efficiently accelerated along the z- 
axis close to the black holcl^^ This would obviously allow acceleration beyond the 
asymptotic flow speed in ideal MHD 7 ~ am (e.g., leiU^. As shown by Beskin & 
Rafikovim for a quasi-monopole magnetic field configuration, this could happen in 
a narrow boundary layer of thickness Ar close to the generalized light surface 

where \E\ — |^^^^ x i?| ~ \B\. Here, 77c > characterizes the strength of the longitu- 
dinal current in terms of the Goldreich- Julian current, jn = (1 — ric)JGJ- ^oi small 
longitudinal currents, i.e. rjc ^>- 1, the light surface coincides with the classical light 
cylinder 

n. ^ rL = ^ . (30) 

Sampling an electric field of strength ^ |-B|, electron acceleration (in the absence of 
losses) is then approximately described by 

^ - ^ (31) 

ar m^c^ 

The maximum electron Lorentz factor that can be achieved in the boundary layer 
of width Ar thus becomes 

7^§^Ar. (32) 

ar 

As the thickness of the layer depends on the charge density, i.e., 

Ar~^, (33) 

where A = n^/nQj ^ 1 is the multiplicity factor and uq] = il_B/(27rce) the 
Goldreich-Julian density,'^^ the maximum Lorentz factor becomes 

7 - ^m (34) 
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where am — B^ / (STmemeC^) is the MicheP^ magnetization parameter. Thus, in a 
narrow boundary layer close to the light cylinder tl almost all of the electromagnetic 
energy may be transformed into the kinetic energy of particles. 

3.3. Centrifugal acceleration close to the light cylinder 

A somewhat different approach, that does not rely on direct electric field acceler- 
ation but takes inertial (centrifugal) and radiation reaction effects explicitly into 
account, has been explored in the literature as welll"^l"^l"^l"^l"^l^^°l^l^^^l 
This approach assumes that the plasma density is in fact high enough to ensure that 
the parallel electric field component is effectively screened out up to distances (of at 
least) very close to the light cylinder, and that the magnetic field structure is such 
as to allow a significant fraction of the electromagnetic energy to be transformed 
into the kinetic energy of particles through magneto-centrifugal effects close to tl. 
This goes along with the requirement that the longitudinal electric current is small 
enough. Consider for illustration a charged test particle, co-rotating with the mag- 
netic field and gaining energy while moving outwards. Under most circumstances, 
an energetic electron will quickly lose its perpendicular energy component due to 
strong synchrotron losses and may thus be considered as just sliding along the field 
line. The motion of such a particle along a field line may be conveniently analyzed 
in the framework of Hamiltonian dynamics: The Lagrangian L for a particle with 
rest mass toq moving along a (two-dimensional) relativistically rotating field line 
(angular velocity Q = c/rL = const.) is simply given by 

L = -moc^l-v'Jc'-vl/c'r/'. (35) 

where Vr = r and v^ = ilr + rj^ (note, that as the field is considered to be 
swept back, B^/Br < 0). As L is not explicitly time-dependent, the Hamiltonian 
H = fP — L, with P = dL/dr the generalized momentum, is a constant of motion 
(Noether's theorem). Using the above relations, one easily finds (cf. also ref.'^^^l) 

H = ■ymoc^ I 1 ^W0 J = const. , (36) 

where 7 = (1 — vf/c^ — w5/c^)~^'^ is the Lorentz factor of the particle. This can 
be generalized to the case where B and f2 are inclined at an angle 6 by replac- 



ing r — )• (rsinfl) on the rhs of eq. (36) and, correspondingly, in Ferraro's law of 
isorotation. Note that the above expression might have also been obtained from 
the standard Bernoulli equation (energy fiow conservation). Provided the product 
fBfjj/ Br remains sufficiently small, this implies 

Thus, for a particle approaching the light cylinder r — > tl, the Lorentz factor 
may increase dramatically as long as co-rotation holds. However, even in the single 



September 26, 2011 0:20 WSPC/INSTRUCTION FILE review 



Nonthermal Processes in Black-Hole-Jet Magnetospheres 23 

particle approach the Lorentz factor cannot become arbitrarily large, as from a 
formal point of view the validity of the approach demands that at least Wc < i^ (see 
also the breakdown constraint below) ,1123 where uJc — sB/{'-fmoc) is the relativistic 



gyro- frequency. Using Eq. ( 36 1 and the definition of 7, the characteristic acceleration 
time scale can be expressed as^^ 

7 1 

tacc — — — f, ^ -1/4 T72 ' '- •' 

valid for 7 S> 1, where to = l/(7o [1 ~ ^o/'"l]^) depends on the initial injection 
conditions. Under more realistic circumstances, achievable particle Lorentz factors 
are expected to be limited either (i) by radiative losses (curvature or Compton), 
(ii) by the breakdown of the (single particle) bead-on the wire approximation (rel- 
evant for, e.g., protons), or (iii) - if we consider an ensemble of particles - by the 
particles' inertia overcoming the tension in the field lines (i.e., breakdown of plasma 
CO- rotation, e.g., see ref .1^211) . In the simplest case the breakdown (ii) corresponds 
to the situation where the Coriolis force becomes comparable to the Lorentz force. 
It is tantamount to requiring that the inverse of the relativistic gyro-frequency Uc 
remains smaller than the acceleration timescale.'i^il This translates into an upper 
limit on achievable Lorentz factors of 

— °'--(^)"'(S)"(isa^)" <-) 

where B(ri^) denotes the field strength at the light cylinder tl. 

3.4. Magnetic reconnection at the jet base 

Efficient magnetic reconnection (annihilation of magnetic fields) could possibly take 
place close to the black hole if there are suitable field regions with opposite polarities. 
In the reconnection process, the energy stored in these fields is released as kinetic 
energy and heat (increase in entropy) (cf. ref.l^^ for review). This requires finite 
resistivity (i.e., breakdown of ideal MHD and the frozen-in condition) and becomes 
important in regions with large magnetic field gradients (for application to AGN 
see, e.g., refsl23I2Zl). By Ampere's law, V x B ~ {in/cjj, these large field gradients 
are associated with large current densities j ("current sheets"). Thus, even for a 
relatively high electrical conductivity Ce, significant ohmic losses (j'^/oe) might 
occur. 

The rate of magnetic reconnection is controlled by the geometry of the dissipation 
region. In the simplest (Sweet-Parker) picture, see Fig. pj the inflow of plasma 
and magnetic field into the reconnection region at speed Vj. is balanced by an 
outflow at speed va- The characteristic outflow speed for a magnetized plasma 
is the Alfven speed va = i3/(47rncTOp)^/^. Mass conservation approximately implies 
L Vr ^ dR Va, i.e., Vr ^ {dR/ L)va- From the (steady state) induction e quation, 
the width can be estimated to be dR ^ X^/vr, where A„j = c^/(47rcre) oc 1/Rm 
is the magnetic diffusivity {Rm ^ "f^AL/Xm 3> 1 the magnetic Reynolds number). 
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measuring the ability of a magnetic field to diffuse through a plasma. This coni- 

1 /2 

bines to Vr ^ va/RtL <C wa, so that the reconnection process would only proceed 
at very slow speeds, the classical problem associated with (Sweet-Parker- type) re- 
connection. 

However, plasma micro-turbulence or wave-particle interactions might lead to an 
enhanced ("anomalous") resistivity (increasing A), allowing faster reconnection rates 
with Vr — >■ va- Such a situation has been envisaged to take place at the interface 
between a black hole magnetosphere and a coronal wind, where poloidal magnetic 
field reversal may occur, see Fig. |g| i28 | i29 | rpj^g maximum rate of magnetic energy 




i«' 




'mmm 



Fig. 8. Left: Schematic drawing of the assumed magnetic field geometry. Acceleration could occur 
in the magn etic reconnection region at the Y-type neutral zone. From de Gouveia Dal Pino and 
LazariarBZll Right: Illustration of a simple reconnection region. Two oppositely directed magnetic 
fields in a plasma are carried towards the neutral line at speed Vr over a characteristic length scale 
L. There is a layer of width dR in which the field reconnects. Reconnected field and plasma are 
expelled at speed va- 



that can be extracted from the reconnection zone in the corona (above and below 
the disk) would then be on the order of 



Pb - (BV87r)w^(47ri?)L - (B'^ /8Tr)vAi'iTrR)dR, 



(40) 



where B is the magnetic field at the reconnection zone, R ^ R^ is of the order of 
the inner radius of the disk, and va ^ c for a highly magnetized plasma. In analogy 
to Ampere's law, we can roughly estimate the width dR of the current sheet for 
which the resistivity must be anomaloufP^ 



,„ cAB m„c^ 

dR 200 -^- , 

ATTUcZevth.c ZeB 



(41) 



where Wth.c 
X-ray {kT 

G(cf. 



i {kT/nipY^^ ^ c/100 is the typical thermal velocity of the ions in the 
100 keV) emitting corona, and where AS - 2B - 10^ (IO^Mq/M)^/'^ 



2.4) denotes the change of the magnetic field across the reconnection region. 
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This gives a relatively modest power of 

/ M \ 1/2 
Pb - 5 X 10^\j^^^) erg/s. (42) 

Yet, while not sufficient to account for radio-loud AGNs, this power could perhaps 
be sufficient to explain the observed radio luminosities of radio-quiet AGNs (see 
Fig. 3 in refp2l). 

During the reconnection process, efficient acceleration of supra-thermal particles in 
the inductive electric field E^. = Vr x B/c perpendicular to the (two-dimensional) 
reconnection fields can occur as long as their Larmor radii do not exceed the width 
of the current sheet P^ Ignoring energy losses, this would suggest that in the con- 



sidered scenario (cf. eq. 41 ) relativistic energies up to 7p ~ 2 x 10^ (protons) and 
7e ~ 4 X lO'"^ (electrons) might become achievable (cf. also ref.li^SI for possibility of 
further acceleration). 

The details of the motion of charged particles in current sheets are non-trivial,12ni 
and obtained results are dependent on the (time-dependent) sheet structure and, 
e.g., the strength of the linking magnetic field perpendicular to the sheet. Pos- 
sible particle spectra discussed in the literature for fast magnetic reconnection in- 
clude quasi-monoenergetic distributions (e.g., ref.l^^ and power-law particle spectra 
n(7) ex 7~* with s ~ 1 (e.g., ref.^^^) or s ~ 1.5 {lei^^^. 

4. RADIATIVE PROCESSES 

The interactions of energetic particles with (photon or magnetic) fields and/or ambi- 
ent matter usually provide important efficiency constraints for the above mentioned 
acceleration mechanisms. Some of these loss processes and their possible relevance 
are discussed below. 

4.1. Synchrotron Radiation 

A charged particle with Lorentz factor 7 and charge number Z spiraling in a mag- 
netic field B emits synchrotron radiation with a mean frequency 

{ly) = 0.31 ly,, (43) 

where Vcs = (3/47r)r2e7^ sin0, flf, = ZeB / (rriQc) is the gyro-frequency and 9 the 
pitch angle between the magnetic field and the velocity vector. The energy loss 
rate, —dE/dt, or the power emitted by a single particle is 

so that the characteristic cooling timescale tg = E/\dE/dt\ = ^rriQC^ / P^s becomes 
— '"^^'^ "■520r!^Vri^V . \, sec. (45) 



2 Z^e^-iB^ain-'e \m^ J \ B J Z^-ism'-Q 
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In terms of the mean frequency (v) , one obtains 

t, ~ 5.9 X 10"Z-'^/2(mo/me)^/2(Ssin6l)-^/2 (i,)-i/2 ggp_ (45) 

The characteristic cooling time of a proton, for example, emitting synchrotron radia- 
tion at (v) = 2.4x10^6 jjz (energy of 1 TeV) is thus of order t^ ~ 64 (B siii9)-^/'^d. 
The synchrotron spectrum (power per unit frequency) I{i') emitted by a single 
particle is a power law I{iy) ex ly^'^ for v ^ z^cs, and decreases exponentially 
for v ^ Vc- For most practical purposes, it can be reasonably approximated by 
I{v) ex {v/vcs)°'^e-Kp{-v/vcs)- 



Equation (451 shows that for almost all pitch angles 9, energetic electrons are ex- 
pected to quickly radiate away their perpendicular momentum components in the 
strong magnetic fields close to the black hole. Electrons in the black hole magneto- 
sphere may thus be regarded as mainly occupying the ground (n — 0) Landau state, 
moving quasi-one-dimensionally along the field lines. In order to emit synchrotron 
radiation, particles in the ground Landau state would have to be excited to higher 
Landau levels (by acquiring non- vanishing pitch angles). This is not excluded but 
may become possible through, e.g., Compton scattering or non-resonant diffusion 
in pitch anglesEMEMHH 

4.2. Curvature Radiation 

In the black hole magnetosphere, and in particular closer to the light surface, the 
curvature of the magnetic field may no longer be negligible, so that particles moving 
along the fields might efficiently lose energy due to curvature radiation. In analogy 
to synchrotron radiation, curvature radiation can be described as emission from 
relativistic charged particles moving around the arc of a circle chosen such that the 
actual acceleration corresponds to the centripetal one, i.e. with curvature radius 
Re = "/moc^/{ZeBsm6) (e.g., ref.^^^. The critical frequency where most of the 
radiation is emitted is given by 

-c ^ ^1' , (47) 

which for a curvature radius of, e.g. Re = rs — 3 x lO"'^'^(Af/lO^M0) cm yields 
i/c =i 2 X 10^^ (7/101")^ Hz or a curvature photon energy of about 1 {-f/10'^°f TeV. 
Like synchrotron emission, the spectrum /(t^) produced by curvature radiation of a 
single particle is a power law I^v) ex v^^'^ for i> <$:. i>c, and decreases exponentially 
for V ^ Vc- The energy loss rate or total power radiated away by a single particle 
(charge number Z) is 

The characteristic cooling timcscalc t^ = ^mQ(? j P^ thus becomes 

^.^180|("^)1. (49) 

Z^ \vne} T 
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In the absence of other damping mechanisms, a balance with coohng dE/dt = Pc 
(using eq. |27| ) would imply that direct electric field acceleration may account for 
particle Lorentz factors up to (cf. refP^^ 



provided the potential (gap height h, spin a, cf. eqs. 25 26 ) is such that these 
Lorentz factors can in principle be achieved. Note that eq. (50 1 does not depend on 
the particle (e.g., electron or proton) mass. The maximum energy of the emitted 
curvature photons becomes 



3 / n.r \ 1/2 / o X 3/4 . p . 1/2 . , . 3/4 



3 c/it' max _ 2 a3/4 ( ^^ V ( ^^ \ 



2 i?c ^ "" \IQ^Mq) IvlO^Cy \rH) Vh, 

(51) 



RcV'U h 



TeV. 



Equation ( 50 ) suggests that the maximum energy for a charged particle of mass 
jjinip may exceed ~ 10^^/i eV only in massive black hole environments. The magnetic 



field strength typically scales as B ex M~i/^ (cf. § 2.4), so that efficient UHE proton 
acceleration would require a very massive source. Equation (51 ) then suggests that 
this should be accompanied by significant VHE 7-ray emission. 
As an illustration, let us suppose that particles are efficiently accelerated in the 
magnetosphere of a massive black hole, with their maximum energies only limited 
by curvature losses. As noted above, this could lead to detectable emission in the 
TeV energy band. For the gap to exist, the charge density has to be smaller than 
the Goldreich- Julian value, cf. eq. ([23]), i.e., n < iiqj = flB/{2Trec) (provided this 



density can be supported by the fields). Assuming a quasi mono-energetic particle 



distribution, the total curvature output at energy ec, eq. (51 1, is then of order 
Lc — n AV Pc- If acceleration takes place in a gap (of thickness h, with associated 
volume element AV ~ nrjr'^h, where 77 < 1 is a geometrical factor), we have at 



maximum P^ = dE/dt = e^ec/h (cf. eq. 27 ). This would imply a maximum VHE 
luminosity of the order of {e.g.,^ 






n (nrjrlh) ( ^ ) (52) 



8 X W''?] 



2/ ,^ \2/,s2 



fj^\ fj^y ( M y /^y erg 



It is interesting to note that the ratio Lc/P^z, where Pbz denotes the maximum 



Blandford-Znajek power, eq. (28), scales as (h/rg)^ (see also ref.^). 



4.3. Inverse Compton Scattering 

Compton scattering of photons by relativistic electrons (Lorentz factor 7 3> 1) is 
usually called inverse Compton (IC) scattering. The IC power for single scattering of 
an isotropic photon distribution with energy density f/ph ~ /i J " di/ i^ripij(i^) 
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[erg/cm'^] in the Thomson Umit hv ^ meC^ /j (i.e., for neghgible energy transfer in 
the electron rest frame) is 

Pic^~ = l'JTCi'U^^, (53) 

which gives a characteristic electron cooling timescale of 

tic - 3.1 X 10^—— sec . (54) 

lUph 

The photon scattering rate, i.e., the total number of incident and up-scattered 
photons per unit time is cuTn^h, where riph = Upi^/h {v). From P^ = cuttIpyJi (vi) 
one deduces that the mean scattered photon frequency is 

('^i) = ll' (^) (55) 

where (i^) is the average incident photon frequency. From the kinematics, the maxi- 
mum (peak) photon energy after scattering is vi ~ A^y'^v for radiation with incident 
frequency v. Taking the Thomson condition into account, we have hvi ~ jnieC^- If 
Compton recoil cannot be neglected (i.e., if the Thomson regime does not apply), the 
full Klein-Nishina cross-section must be used, resulting in a reduced power output. 
Based on QED considerations, the total cross-section is given by the Klein-Nishina 
(KN) formula, for which two useful limits are 



(Jt for hi^r <C rrieC 
3 f rn^c 
8 \ hi^r 



, 2hvr \ 1 
In 



, , for hvr ^ rrieC . (56) 

^nieC^ J 2 J 

where Vr is measured in the rest frame of the electron. In the extreme KN regime, 
the energy loss rate for a single electron is found to b^^^Si 

/A-yhv 11' 



dE 3 2 51,-1 f J nph{iy) 



In 



\7TT,eC^ 6 



(57) 



and is thus only weakly (logarithmically) dependent on the electron Lorentz factor 
7. However, in contrast to the Thomson case, the photon now carries away a sizable 
fraction of the electron energy, so that the electron loses its energy in discrete steps. 
The spectrum emerging from single inverse Compton scattering generally depends 
on both, the seed photon spectrum and the relativistic electron distribution. If the 
range of the seed photon spectrum is very narrow compared to the electron dis- 
tribution (assumed to obey a power-law rig (7) oc 7^* between 7i„in and 7max), for 
example, the output spectrum is a power-law v~y^^^il'^ over a dynamical range of 
(7max/7min)^, with indcx determined by the electron power-law index. On the other 
hand, if the electron distribution is very narrow compared to the photon distribu- 
tion, the output spectrum resembles the seed photon spectrum shifted in energy by 
a factor ^ (7^)- 

If electrons are rotationally accelerated in the magnetosphere of a black hole embed- 
ded in an ambient photon field, Compton losses may significantly reduce achievable 
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particle energies. In the Thomson regime, for example, inverse Compton upscat- 
tering of accretion disk photons will lead to particle energy losses on a timescale 



he c>c 7~^ that decreases faster than iacc oc 7"^' ^, eq. (38 1, and so introduces a nat- 
ural limitation. Balancing acceleration by cooling, eq. (54 1, implies a characteristic 
upper limit of (e.g., ref.^^^) 



,6 v^ 



3 X 10 772" ' (58) 



assuming co-rotation to hold for such a range of Lorentz factors. Hence, only for 
highly underluminous AGN sources with target photon fields Lt <Si L^dd (Uph ^ 
Lt/[4:TTr1c]), will centrifugal acceleration allow to accelerate electrons to Lorentz 
factors 7 » ioo)ii8|ii9| 

4.4. Photo-meson production 

Photo-meson production can become an important channel by which the ki- 
netic energy of protons is transformed into high energy gamma-rays, electrons 
and neutrinos.'i^ The kinematic threshold for single pion production (p+7 — )• 
p -\- tt'', n -\- TT"*") is given by 

2Epe-y{l - cos 61) = (2m^mp + m\)c^ ~ 2.8 x 10^^ (eV)^ (59) 

where E^ is the proton energy (assumed to be relativistic), e^ = hv is the target 
photon energy and m.,^ the pion rest mass. Thus, in the rest frame of the proton, the 
photon energy has to exceed e' = TO7rC^(H-m7r/2m.p) ~ 140 MeV for the process to 
become possible. Close to the energy threshold, the process proceeds through single- 
pion production, while at higher energies multi-pion production channels start to 
dominate. The cross-sections for these inelastic processes are known from particle 
acceleration experiments. In first-order approximation, they can be expressed as a 
sum of two step-functions with a^ = 0.34 mbarn (= 3.4 x 10^^* cm^) in the interval 
200 MeV < e'^ < 500 MeV and cts = 0.12 mbarn for e' > 500 MeV, and related 
inelasticities K^x — 0-2 and Kp2 — 0.6, respectively.ESH Por the characteristic 
attenuation length we have A ^ l/{n^aKp), where n^ is the photon number density 
[particles/cm'^] , so that the characteristic proton cooling time t^ — X/c becomes 

. / e^ \ /I042erg/s\ f R ^^ 



where rij ~ L/(47r i?^ce^) has been expressed in terms of the source luminosity L. If 
the inner disk surrounding the black hole would be of the standard (geometrically- 
thin, optically-thick) type, then the thermal disk radiation field would be dominated 



by emission from regions close to r^, with energy peaking at e^ ^ 50 eV, cf. eq. (151, 
and thereby provide a suitable ambient target field for protons with Lorentz factors 
above jpt ^ 2 x 10® (50 eV/e-y). Neglecting curvature losses and field screening, gap- 



type particle acceleration of protons up to ultra-high energies ~ 10 eV, cf. eq. (261, 
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would then only be possible if mass accretion occurs at low rates corresponding to 
disk luminosities Ld ~ O.OlLEdd (cf. refP2l). 



4.5. Bethe-Heitler (proton-photon) pair production 

At energies below the threshold for photo-meson production, the main channel for 
inelastic interactions of energetic protons with ambient photons is the direct produc- 
tion of electron-positron pairs, i.e. the process p+7 — > e++e^+p, e.g.refsEMSIlThe 
minimum energy (threshold condition) required for this process to become possible 
is 

7„e..(l - ^p cos 61) > 2 (nieC^ + —m,^c^) , (61) 

nip 

where 7p = Ep/nipC^ is the proton Lorentz factor and e-y the ambient photon energy. 
Thus, for head-on collisions the process is energetically allowed when 7pe^ > rrieC^. 
The maximum energy of the resultant electron (positron) is determined by the kine- 
matics and (for nieC^ <C ^p£-y ^ nipC^) given by (cf. inverse Compton scattering) 
-Ee,max = 47^6/. The cross-section for this process (often referred to as Bethe-Heitler 
cross-section) is about two orders of magnitude higher than the one for photo- meson 
production, but on average only a small fraction of the proton energy (~ me/mp) 
is lost. This is different to photo-meson production, where a proton transfers on av- 
erage '^ 13% and more of its energy to the secondary products. As a consequence, 
photo-meson production usually becomes the dominant energy loss channel above 



the threshold eq. ( 59 1 



A variant of the Bethe-Heitler pair production process in the ergosphere of a maxi- 
mally rotating, supermassive black hole (referred to as Penrose Pair Production) has 
been suggested already some 30 years ago as an efficient way for energy extraction 
from rotating Kerr black holestSMSSI Upon entering the ergosphere, photons with 
energies of some tens of MeV (as might be produced in a hot ADAF inner disk via 
decay of neutral pions from pp-coUisions, ref.'^^Sl) could get boosted (blue-shifted) 
by a factor of up to ^ 30 (for extreme Kerr black holes) and so approach GeV 
energies. These photons could then interact with the (low-energy) protons moving 
on marginally stable orbits deep within the ergosphere to produce electron-positron 
pairs. From the threshold condition, minimum photon energies of only 2mcC^ would 
be required. However, in order to become a Penrose Process, the recoiling protons 
must be injected on negative energy orbits (and pass through the event horizon), 
with the escaping pairs ejected with an energy boost picked up from the rotational 
energy of the Kerr black hole. Kinematically this requires the blue-shifted photon 
to have an energy comparable to the rest mass energy of the proton, giving GeV 
energies to the escaping electron-positron pairs. While this process may in principle 
be possible, detailed Monte Carlo simulations suggest it to be less important as 
most of the pairs apparently do not manage to escape.C^ 
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4.6. Inelastic proton-proton collisions 

Relativistic protons and nuclei can produce high energy gamma-rays in inelastic col- 
Hsions with ambient gasP^ The neutral 7r°-meson (with main ~ 99% decay mode 
into two photons, e.g., p+p — > p+p+n'^, tt*^ — > 27) provides the main channel of con- 
version of kinetic energy to high energy gamma-rays. For Tr'^-production to become 
kinematically possible, the kinetic energy of the bombarding proton must exceed 
Eth = 2mT^c^{l + rnTr/[4:mp]) ~ 280 MeV. The cross-section and the inelasticity for 
pp-coUisions both depend weekly on the energy, i.e., app ~ 3.4 x 10~^^ cm^ and 
total inelasticity /^ ~ 0.5. The characteristic cooling timescale tpp ^ ^/{(^ppnpf^c) 
for a relativistic proton thus becomes 

/ lo^cm"^ \ 
tpp ~ 2 X 10^ ( — ) sec , (62) 



n. 



V 



where Up is the ambient gas density. While tpp is only weakly dependent on the 



energy of the proton, the cooling timescale for p7-reactions, cf. eq. (60), is expected 
to decrease with (increasing) proton energy (i.e., lower photon energy). Under such 
conditions, pp-interactions can dominate the cooling only below a certain proton 
energy (cf. also ref.^^^ for an application to AGN). 

In the simple delta- functional approximation,'I5flithe neutral pion takes away a mean 
fraction /^o ~ 0.17 of the kinetic energy of the proton £"^0 = /TrO-Ep.kin, and this 
energy is equally distributed among its decay products. The shape of the photon 
spectrum is then similar to the shape of the parent proton spectrum, but shifted in 
energy by a factor /^o . 



4.7. jj -Absorption 

Photon-photon interactions can produce electron-positron pairs (7 + 7 — )• e+ -|- e^ ) 
and thereby lead to a suppression of the VHE photon flux from a source, once the 
kinematic threshold condition for pair production {hv){hvt) > 2(TOeC^)^/(l — cosa) 
is satisfied (a being the angle between the incident directions of the photons). 
TeV photons, for example, can thus react (head-on) with infrared target photons 
of energy £< = hvt ^ 0.3 eV. The electron-positron pairs produced in energetic 
collisions will be highly relativistic and tend to move in the direction of the initial 
VHE 7-ray. 

For a fixed photon energy e, the pair-production cross-section a-^-y is a function of 
target energy ct P^ If we denote by s = i?^/(4TOgC'*) the normalized centre-of-mass 
energy of the photons squared {E"^ = 2eet[l — cos a]), then starting from zero at the 
threshold energy s = 1, the cross-section rises steeply to a maximum ~ (i7t/4) at 
s ~ 2 and decreases approximately oc 1/s ex 1/et towards higher energies, i.e., 

a^,{s)^l^[H4s)-l] (63) 
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for s ^ 1. In an isotropic background radiation field, TeV photons thus interact 
most efhciently with infrared target photons of energy 

/lTeV\ 



et-1 { ^— j eV. (64) 



The optical depth, characterizing the absorption of a 7-ray with energy e moving 
through a target photon gas of spectral and spatial distribution n{r, et) in a source 
of size R, can then be calculated from 

r(e) — dr dct n^{r^et) (T-y^(e, €4, a) (1 — cosa). (65) 

Jo -'2(mec2)2/[e(l-cosQ)] 

For a homogeneous source and a quasi-isotropic target photon field, the optical 
depth become^ 

/>oo 

T(e) ~ i? / dtt n^{et) a-^j , (66) 

where n-y{et) is the differential target photon number density. Because the optical 
depth for a gamma-ray of energy e is essentially determined by a relatively narrow 
band of target photons (centered at s = 2), a useful order-of magnitude estimate is 
given by r(e) ~ 2.5{<7t /4)Retn-y{2et)!^^ln terms of the luminosity Lt above energy 
Ct on 1/e, the optical depth may thus be approximated by 

When applied to AGNs, this suggests that only in highly under- luminous {Lf <C 
^Edd) and massive (i? ex rs oc M) sources, VHE gamma-rays may be able to 
escape unabsorbed from the vicinity of the central supermassive black hole. As an 
application, let us suppose that the inner accretion disk is of the high-temperature 



{Te ^ 5 X 10^ K) ADAF type (see § 2.3.4). The dominant part of the radiation 



in the hard X-ray regime {hv ~ kTg ~ 500 keV) would then be due to thermal 
bremsstrahlung. These hard X-ray photons would preferentially interact with target 
photons of energy et ^ 2 MeV (cf. eq. [64]). Although the disk spectrum would be 
suppressed at these energies (as the bremsstrahlung spectrum falls of exponentially 
with exjpl—hiy/kTe]), it may still be strong enough to lead to some non-negligible 
electron-positron pair production. We can roughly estimate the density of pairs so 
provided by balancing the rate of their creation with the rate of their escapejj The 
pair creation rate [particles/sec] is of the order of Ni ~ {4/3)TTR^n^/t\ifo where 
iiife ~ l/(n-y(T-),^c) is the characteristic lifetime of an incident photon. The created 
pairs will escape from the source region at a rate N^ ^ AirR^rLf-c (ignoring Compton 



A full analysis would have to employ the appropriate angle-averaged cross-section a^^, see e.g. 
[T52l 

'I owe this argument to Marek Sikora. 
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scattering). Hence we find rig ^ nia^^R/S. Using n^ ^ Lt/(47ri?^cet), we have 
L^cr^^ ^/ Lt \ VlO"cm\^ /2 MeV^^ 



/ Lt y/lO"cm 



487r2i?3c2g2 \^1040ergy \ R J \ et 



i 77 o / i \ I \ I ivit; V \ - _3 



[particles cm ]. (68) 



Given the conventional AGN parameter space, one may expect the so-estimated pair 
density to be comparable to the Goldreich- Julian density in a fair number of sources, 
implying that a substantial part of the parallel electric field in the magnetosphere 
may be screened. This may not be the case, however, for highly under-luminous 
systems such as M87 and Sgr A*llilllln terms of the accretion rate, the density ratio 
of so-created pairs (in a two-temperature ADAF) to the GJ value approximately 
foUowslS 

Hence, if the accretion rate would become very small, an additional plasma source 
(such as cascade formation in starved magnetospheric regions) would be needed to 
establish a force- free outfiow (where n,, > uq]). 

Note that if a source is sufficiently compact, its VHE 7-ray emission could well 
be shaped by the development of an electromagnetic cascade (photo-production of 
pairs, subsequent inverse Compton scattering etc.).'^^ Also, since both TeV 77- 



absorption and photo-meson production, § 4.4 can interact with similar target pho- 
ton fields, TeV gamma-ray observations could constrain the photo-pion production 
opacity, in particular for sources where rapid VHE gamma-ray variability is ob- 
served. When applied to TeV blazars, this suggest that the expected flux of high 
energy neutrinos from p7-interactions is much less than the observed VHE gamma- 
ray flux, making them less promising targets for neutrino observations.^ 



5. APPLICATIONS 

5.1. Low-luminous galactic nuclei 

According to the considerations above, nearby under-luminous and non-aligned 'ac- 
tive' galactic nuclei emerge as prime candidate sources where non-thermal VHE 
processes, occurring close to the central black hole, may become observable and 
thereby allow to probe the most violent region of the central engine: If the source 
is sufficiently under-luminous, VHE gamma-rays may be able to escape significant 
absorption. The source must then, however, be relatively nearby in order to be 
within the reach of current or upcoming telescope sensitivities. Finally, if the jet is 
sufficiently misaligned, the nearby black hole VHE emission may not be swamped 
by Doppler-boosted jet emission. This makes sources such as Sgr A* or the nearby 
radio-galaxies M87 and Cen A to interesting candidates. 
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5.1.1. The galactic centre source Sgr A 



* 



Given its proximity, the compact radio source Sgr A* in the central region of our 
Galaxy may represent the most promising target for studying particle acceleration 
processes near the event horizon of a supermassive black hole. Believed to host a 
black hole of mass Mbh — (3 — 4) x 10^ Mq^ its extraordinary low bolometric lumi- 
nosity (< 10~*LEdd) is expected to make the innermost region transparent to 7-rays 
with energies up to ~ 10 TeV.'^^ Estimates for the accretion rate in Sgr A* are in 
the range M ~ (0.1 - 4) x IQ-^Mq/ji^^IMMMso that m = M/M^dd ~ 5 x IQ-^. 
This suggests that the magnetic field in the innermost region could possibly be as 
high as S - 5 x 10^ G, cf. p^ 



Early VHE observations of the Galactic Centre region by Imaging Atmospheric 
Cherenkov Telescopes have led to the detection of a steady, point-like source of 
VHE gamma-rays at the gravitational center of our Galaxy, coincident with the 
position of at least three counterparts, the supermassive black hole Sgr A*, the 
supernova remnant (SNR) Sgr A East (at a projected distance of 3.7 pc to Sgr A*), 
and the pulsar wind nebula G359.95-0.04 (at a projected distance of 0.4 pc to 
Sgr A*).J60 I6i|i62|i63|i64j Recent H.E.S.S. resuhs based on high precision pointing, 
however, now seem to exclude Sgr A East as the dominant source of the observed 
VHE gamma-rays .li^ On the other hand, both G359.95-0.04 and Sgr A* remain 
promising VHE counterpart candidates, as models exist that could explain the 
observed emission extending beyond 10 TeV (e.g., see reiU^ for G359.95-0.04). 
In the case of Sgr A*, for example, the origin of the observed TeV emission has 
been related to (i) non-thermal processes in the black hole magnetosphere itself,'^^ 
(ii) pp-interactions of escaping high energy protons in the surrounding dense gas 
environment jiMIMilll or possibly (iii) electron acceleration in the termination 
shock of a wind emerging from the innermost parts of the disk.l^ 
Suppose that particles could be efficiently accelerated in voltage gaps close to the 



black hole horizon. Then assuming a rapidly spinning black hole, eqs. (25 1 and 
(50) suggest that protons may possibly be able to reach energies up to ~ lO^'^ eV 



(7p ~ 10^).'i^lf so, then curvature emission is expected to peak around ~ 10^° eV, 



eq. (47), i.e., well below the TeV domain. However, detectable TeV 7-ray emission 
may still be possible due to other radiative channels, and at least three different 
scenarios have been proposedP^ 

(1) Photo-meson interactions: If protons are able to reach energies of ^ 10^^ eV, 
they can start to interact with infrared soft photons, cf. eq. ( [59| . Although the 
mid-to-near infrared source in Sgr A* is faint with Lir ~ 10^^ erg/s, it appears 
sufficiently compact rm ~ 10^"^ cm (e.g., refsIii3iZ2EZ3) to ensure a high enough 
seed photon density n^ ~ (10^^ — 10^"^) cm~'^ for effective interactions with pro- 
tons. The mean free path of protons through such a photon field, however, would 
be comparatively large, i.e., A ^ l/{apjnjKp) ^ (10^^ — 10^^) cm, cf. [4.4 This 



suggests that only a modest fraction rip/A < 0.01 of the energy in protons can be 
converted into secondary particles. Thus, in order to account for a VHE flux level 
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of Lj ~ 10'^^ erg/sf^^ an injection power up to ^ 10'^^ erg/s in high energy protons 
would be required. 

(2) Proton-proton interactions: If protons are not accelerated to energies ^ 10^* eV 
in Sgr A* (as e.g. expected in the case of centrifugal acceleration), the efficiency of 
photo-meson interactions would be significantly reduced. Interactions of energetic 
protons with the ambient plasma could then become the main channel of VHE 
7-ray production. The associated cooling time scale for pp-interactions is of order 



tpp ^ 2 X 10^ (10^cm~3/np) sec, cf. eq. (62 1, for typical densities n^ -- (10^ - 10®) 



cm~^ in the inner regions of the accretion fiow near Sgr A*, cf. eq (24 1. If the ef- 
ficiency of gamma-ray production is determined by the ratio e^ of accretion time 
scale ta ~ r/vr ~ 10^ sec (with Vr ~ avf, vj the free fall velocity) to cooling 
time scale jli^ then the efficiency of converting the energy of accelerated protons 
to secondaries may be as small as e^ ^ 10~^. Hence, a high injection power of 
L^/cr ^ (lO'^^ — 10^°) erg/s would be required, very close to what may actually be 
provided by accretion. 

(3) Inverse Compton scenario: Upscattering of soft photons by accelerated elec- 
trons may provide the most economic way to produce VHE 7-rays. If the infrared 
source represents the main target field, inverse Compton upscattering will proceed 
in the Klein-Nishina regime, requiring electrons energies ~ 10 TeV to account for 
the observed TeV emission. In order to avoid suppression by dominant electron 
synchrotron losses, particle acceleration must then essentially proceed along the 
ordered magnetic field. Both, centrifugal and/or gap-type particle acceleration are 
potential candidates to this end. If gap-type particle acceleration would be limited 
by curvature losses, for example, up-scattering could result in VHE 7-rays with en- 
ergies Ej ^ 'je'iTiec'^ ~ (10^'* — 10^^) eV, cf. eq. (50). These energetic photons will, 



however, not be able to escape the source, but instead interact with far-infrared pho- 
tons to produce electron-positron pairs. Synchrotron radiation and inverse Compton 
scattering by these secondary electrons will lead to a re-distribution of the initial 
gamma-ray spectrum. Figure [9] shows a model calculation of the expected spectral 
energy distribution within such a scenario. 

5.1.2. The radio-galaxy M87 

Located at a distance of d '-^ 16 Mpc (redshift of z = 0.0043), the giant ellip- 
tical Virgo-cluster galaxy M87 represent another promising candidate for magne- 
tospheric VHE 7-ray emission. M87 hosts one of the most massive black holes 
Mbh — (3 — 6) X 10^ Mq in the Universsii3iI3 and shows a prominent, non-aligned 
jet detectable from radio to X-ray wavelengths. Despite its huge black hole mass, 
M87 is not a powerful source of radiation: Based on its estimated total nuclear 
(disk and jet) bolometric luminosity of Lboi ^ 10^^ crg/s or less,'^^ M87 is highly 
underhiminous with /g < 3 x 10^®, where l^ = ^boiZ-^Edd £^nd L^dd is the Edding- 
ton luminosity. This has led to the proposal that M87 is a prototype galaxy, where 
accretion occurs in an advective-dominated (ADAF) mode characterized by an in- 
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Fig. 9. Possible broadband SED for Sgr A* produced by electron curvature (1st peak) and inverse 
Compton (2nd peak) radiation. The thin solid line represents the gamma-ray spectrum due to the 
accelerated (primary) electrons, the thick solid line shows the redistribution of the ene rgy after the 
passage through the infrared source. Figure adapted from Aharonian & NeronovED (reproduced 
by permission of the AAS). 



trinsically low radiative efRciencyliMIZlEini HST observations of M87 have revealed 
superluminal motion of jet components at ~ 0.5 kpc from the central black hole, 
indicative of bulk flow Lorentz factors r^ ~ 6 and a jet orientation oi 9 ^ 19° to 
the line of sight (ref.'I^ see also ref.^^ for larger 6 based on 43 GHz radio ob- 
servations), suggesting that M87 is a non-blazar jet source, characterized by only 
moderate Doppler factors D = l/[rf,(l — (3cos9)] ~ 3. 

Despite this, recent TeV observationsIEEl have demonstrated that the 7-ray spec- 
trum of M87 extends beyond 10 TeV and is consistent with a relatively hard power- 
law (with spectral index a ^ 1.2, where Si, ex i^~"), see Fig. 10 The TeV output is 
relatively moderate, with an isotropic TeV luminosity of some lO^'' erg/s. Significant 
variability (flux doubling) on time scales of Atobs — (1 — 2) days has been found, the 
fastest variability observed in any waveband from M87 so far. The observed spectral 
and temporary VHE 7-ray characteristics of M87 have been proven difficult to ac- 
count for within classical jet models (cf. ref.l^^for review) and renewed the interest 
into magnetospheric emission models EEEHMI This interest has been strengthened 
by the recently detected correlation between radio VLBA (probing scales down to 
some tens of r^) and TeV 7-ray data during a flaring state of M87 in 2008, see 
also Fig. [ipl There are also indications that a simple extrapolation of the Fermi 
(200 MeV-30 GeV) high energy spectrum under-produces the TeV flaring state,^^ 
supporting the possible appearance of a new (variable) component at the highest 
energies. 

If the observed TeV emission would indeed be due to magnetospheric processes, one 
may expect the light travel time across the source of ~ 0.2 days to provide a lower 
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Energy (TeV) 



Fig. 10. The differential very high energy (VHE) spectra of M87 as obtained by H.E.S.S. during 
a low (in 2004) and a bright (in 2005) state, ranging from ~ 400 GeV to ~ 10 TeV. The indicated 
fits give power law photon indices of (o + 1) = 2.62 ± 035 (2004 data) and (a + 1) = 2.22 ± 0.15 
(2005 data). No variation in spectral shape is found within errors. For comparison: One Crab unit 
at 1 TcV corresponds to a flux of ~ 3 X 10^"'^"'^ s^-"^ ci 



Figure from Aharonian et alP 



limit on the possible variability time scale at VHE energies. This is obviously close 
to what can be probed with existing Cherenkov arrays, and is certainly within the 
reach of future instruments like CTA. The observed day-scale variability already 
implies, e.g., that the (generalized) hght surface, eq. (29 1, must be very close to 



the light cylinder scale of a few Vg, which in the direct acceleration model § 3.2 



equivalent to the requirement of a small longitudinal current. The observed, fast 
variability also excludes, e.g., pp-interactions as dominant channel for VHE 7-ray 
production in M87. 

If the accretion rate in M87 is indeed of the order of the Bondi accretion rate 
M ^ 10~^MEddr^ then possible magnetic field strengths in the vicinity of the 
black hole could be as high a,s B ^ 1000 G, see § |2.4[ This would be consistent 
(if the black hole is sufficiently spinning) with poloidal magnetic field strengths re- 
quired by a Blandford-Znajek-type process (eq. [28]) to account for the estimated 
jet kinetic power of a few 10*'* erg/sEMEZZl 

If electrons are centrifugally accelerated in an ambient photon field of order Lboi, 
achievable particle energies may be as high as 7e ~ (10^ — 10®), cf. eq. (58 1. Similar 



(and most likely even higher) energies might be achievable if efficient gap-type par- 
ticle acceleration would indeed take place, see eq. (27).'3lf the inner disk is of the 



ADAF-type, this could facilitate upscattering (Thomson regime) of (comptonized) 
sub-mm ADAF disk photons, i.e., those soft photons above the synchrotron peak 
i^p ^ 10^* Hz, cf. eq. (20). Upscattering would results in emission at ^ 5 (7e/10^) 



TeV energies and can produce a hard photon spectrum close to what has actually 
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been observed P^ If, on the other hand, the shape of ambient disk photon field 
would be more of the standard disk-type, then secondary inverse Compton cascade 
emission, initiated by internal absorption, could become influential at TeV energies 
as wellP 

An important question for magnetospheric scenarios is whether TeV photons can 
escape from the vicinity of the central black hole. In principle, TeV photons of 
energy e will interact most efficiently with target photons in the infrared regime 



ct — (ITeV/e) eV, eq. (64|. The relevant infrared (IR) photon field needs thus 
to be sufficiently diluted in order to avoid significant pair-absorption. This could 
well happen if the IR emission is dominated by, e.g., emission from the jet (syn- 
chrotron emission) and/or a dusty torus on larger scales. For M87, the observed 
nuclear mid-infrared luminosity has been estimated to be Ljr ~ 10*^ erg/s.^SS The 
corresponding optical depth for TeV photons would thus be r ~ 3 (10^^cm/i?iR), 



eq. (67), assuming the infrared emission to be concentrated in a homogeneous source 
of size i?iR. Hence, for _Rir ~ some tens of r^, TeV photons would be able to escape 
the source. This latter condition is probably difficult to satisfy for a standard-type 
disk where the emission from the innermost part is expected to peak close to the 
infrared regime, cf. § |2.3.2| but it seems well possible for a radiatively inefficient 
accretion flow. If one assumes, for example, that all of the observed nuclear radio 
to X-ray flux in M87 arises in an ADAF (which may appear over-restrictive given 
its jet), then ~ 10 TeV photons would have to be produced on scales ~ Sr^ (for 
a Kerr black hole) and ISrg (for a non-rotating black hole) in order to be able to 
escape the source.'lSZi Recent findings, according to which the infrared emission is 
consistent with optically-thin synchrotron emissioUiliSSI suggests that the real situ- 
ation could be even more relaxed. Hence, there are good reasons to believe that, 
at least for M87, TeV 7-rays may be able to escape from the vicinity of its super- 
massive black hole.'^ This in turn suggests that future high-sensitivity instruments 
like the European CTA project could play a particularly important role in testing 
magnetospheric VHE 7-ray emission scenarios in M87. 



5.1.3. The radio-galaxy Cen A 

As the closest active galaxy (distance d ~ 3.4 Mpc, ref.'^^, the FR I source Cen- 
taurus A (Cen A, NGC 5128) seems another promising candidate worth exploring. 
Being a spectacular example of a gas-rich disk galaxy consumed in a merger with 
a giant elliptical galaxy, Cen A belongs to the best studied extragalactic objects 
(see, e.g., ref.'^^). Radio observations reveal a pecuhar morphology with a sub-pc- 
scale jet and counter-jet, a one-sided kpc-jet, two radio lobes and extended diffusive 
emission. VLBI observations suggest that Cen A is a non-blazar source, its jet being 
inclined at a rather large viewing angle i ~ 50° and characterized by moderate bulk 
fiow speeds Uj ^ 0.5 c.^^^ ■^^'^ Its central black hole mass has been inferred to be in 
the range A/bh = (0.5 - 3) x IO^^MqIIMIEII 
With an estimated bolometric luminosity output of the order of Lf, ^ 10'*'^ erg/s,!!^ 
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Cen A is rather under- luminous (although not highly) and believed to be accreting 
at sub-Eddington rates M ~ lO^'^AfEdd- The apparent lack of a big blue bump 
UV feature seems to indicate that its inner disk is not of the standard-type.l^Sl In 
fact, the disk contribution seems consistent with a hybrid disk configuration where a 
standard disk is truncated at rt and replaced by an ADAF in the inner regions close 
to the central black hole (cf. refsEHlElIElIl). If the inner disk in Cen A would remain 
cooling-dominated (standard disk), magnetic field strengths close to the black hole 



of order i? ~ 2 x 10^ G might be expected, cf. eq. (21). If the disk switches to 
a radiatively inefhcient mode, the characteristic magnetic field strengths may be 
higher, possibly approaching B ^ 10"* G if its black hole mass is at the lower mass 



end, see eq. (22 1 



The possible association of some of the Pierre Auger Observatory (PAO) measured 
UHECR events above 57 EeV with Cen A has recently triggered a number of stud- 
ies analyzing the efficiency of cosmic-ray acceleration in Cen Al ^^^ l ^°° l ^°^ l ^°^ l Based 
on data up to August 2007, the PAO Collaboration initially reported evidence for 
an anisotropy at the 99% confidence level in the arrival directions of cosmic rays 
with energies above ~ 6 x 10^^ gyisosj ^pj^g anisotropy was measured by the fraction 
of arrival directions that were less than ~ 3° from the positions of nearby AGN 
(within 75 Mpc) from the VCV catalog. While that correlation has now decreased 
using the newly available (twice as large) data set, the updated analysis still sug- 
gests that a region of the sky around the position of Cen A has the largest excess 
of arrival directions relative to isotropic expectations .i^^ij jf efficient UHE cosmic- 
ray acceleration would indeed take place in Cen A, this would seem to require the 
operation of an additional acceleration mechanism beyond gap-type particle accel- 
eration, at least in the case of a proton-dominated composition.'222l For even in the 
ideal case where (i) the ordered poloidal field is assumed to be of the order B ~ 10^ 
G, (ii) the black hole to be rapidly spinning a ^ 1, (iii) almost the full induced 
electric potential to be available for particle acceleration, and (iv) radiative con- 
straints to be negligible, achievable cosmic ray energies are not expected to exceed 



^ ~ 2 X W^Z a (M/lO8M0)i/2(M/io-3MEdd)^/^ eV by much, cf. eq. ^. Drop- 
ping (iv) by taking curvature losses into account would reduce achievable proton 
energies to ^ 10^^ eV, see eq. (50). If, on the other hand, a centrifugal-type acceler- 



ation mechanism would be operative, achievable energies may be even smaller, see 
eq. (39). Acceleration of protons to energies beyond a few times 10^^ eV in regular 



magnetic fields close to the black hole would thus seem to be disfavored, even if 
Cen A was in the past in more active stage. The situation is much more relaxed 
in the case of heavier elements like iron nuclei, and this may perhaps be consistent 
with recent PAO indications for an increase of the average mass composition with 
rising energies uptoi?~10^®'^ gyiaosj 

At 7-ray energies, Cen A is the only AGN of the non-blazar type detected at MeV 
(COMPTEL) and GeV (EGRET) energies. Its nuclear spectral energy distribution 
(SED), based on non-simultaneous data, appears to be composed of two peaks. 
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one reaching its maximum at several times 10^'^ Hz and one peaking around 0.1 
]y[gY 206 198 r^j^^ SED bclow 1 GcV has been successfuUy modeled within a simple 
jet synchrotron self-Compton (SSC) framework, assuming Cen A to be a misaligned 
BL Lac object and its first SED peak to be due to synchrotron emission (refl^^S for 
an alternative interpretation, cf. refpSZi). The detection of faint VHE 7-ray emission 
(integral flux of ~ 1% of the Crab Nebula above 250 GeV) up to ~^ 5 TeV has been 
recently reported by the H.E.S.S. coUaborationPSS The VHE spectrum (from ~ 250 
GeV to 5 TeV) can be described by a power law with spectral index a ~ 1.7 ± 0.5. 
No significant variability has been found in the data set. At lower 7-ray energies 
(100 MeV to some tens of GeV), Fermi has recently also detected both the core 
and the giant radio lobes of Cen AEM2121 As it appears, a single population of high 
energy particles seems unable to account for the core flux observed by Fermi and 
H.E.S.S. In fact, a simple extrapolation of the Fermi (power law) spectra would 
tend to under-predict the flux at TeV energies EMHSl This could possibly indicate 
an additional contribution to the VHE domain beyond the conventional SSC jet 
emission, emerging at the highest energies. Obvious candidates within this context 
could be: (i) electron inverse-Compton (IC) emission powered by magnetospheric 
processes (centrifugal acceleration of electrons) in an ADAF environment p^ (ii) 
IC pair cascade emission in the radiation field of a standard disk,E12l provided the 
magnetic field is weak enough for synchrotron cooling to be suppressed (cf. also 
ref .1^^ , and (iii) hadronic (p7) interactions in a standard- type disk environment, 
provided protons can achieve ~ 10^^ eV^^ These approaches are complementary 
in that the leptonic scenario (i) requires a rather low ambient IR photon field close 
to the central source, whereas the leptonic model (ii) or the hadronic model (iii) are 
dependent on rather large UV and/or IR background fields. An exemplary output 
of model (iii) is shown Fig. 



11 



Let us suppose that TeV 7-rays could be indeed produced relatively close to the 
central black hole in Cen A. Then given its putatively small source scale and rela- 
tively high nuclear mid-infrared brightness (isotropic ~ 6 x lO''^ erg/s at hv ~ 0.15 
eV, with evidence for an exponential cut-off towards higher frequencies, see, e.g., 
refsEMIMl), efficient escape of these VHE photons might appear much more chal- 



lenging for Cen A than for M87, cf. eq. ( 67 1 . However, most of the nuclear IR emission 



may well arise on larger scales (note that the mid-IR emission is unresolved on scales 
^ 0.2 pc) and thereby allow sufficient dilution. In fact, the observed nuclear mid-IR 
emission is commonly believed not to be produced closed to the black hole, but 
to be dominated by a non-thermal (non-isotropic!) synchrotron jet component at a 
distance of ~ 0.03 pc (cf. refsEMHIIlIII) and/or a (quasi-isotropic) dusty torus on 
scales ^ 0.1 pc or larger (refsEMIS). If this would be indeed the case, then photons 
with energies of a few to several TeV might be able to escape unabsorbed, while 
absorption features should become apparent towards higher energies. This could be 
tested with future high-sensitivity VHE 7-ray instruments and help to disentangle 
the site of the VHE 7-ray production and the origin of the observed nuclear mid-IR 
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Fig. 11. Photon fluxes for Cen A calculated within a hadronic scenario, assuming p7-interactions 
within a strong UV background photon field. Flux suppression by optical depth effects due to the 
ambient IR field becomes apparent towards higher energies. The different lines show the output 
for different maximum proton energies -Emax assuming either a diffuse or compact IR source, i.e. 
a spatial IR source size of iJjR = 1 pc and emission cut-off in far IR [left], or a source size 0.01 pc 
and no emission cut-off [right]. In order not to exceed the observed VHE fiuxes, proton energies 
~ 10^^ eV seem to be required. Figure adapted from KachelrieB et al 



emission. 

5.2. Passive black holes as TeV j-ray and UHE cosmic-ray sources 

If supermassive black holes would indeed be efficient UHE cosmic-ray accelerators, 
one may expect this to be accompanied by a significant TeV 7-ray output due 
to curvature radiation, e.g. ref.^^^. Observational constraints on the latter may 
then allow to impose physical constraints on the first. Perhaps the most obvious 
candidates for efficient gap-type particle acceleration are massive, but very weakly- 
emitting (dormant or "non-active") sources: 



5.2.1. The Fornax cluster galaxy NGC 1399 

Consider, for illustration, the giant elliptical Fornax cluster galaxy NGC 1399 at a 
distance of ^ 20 IVIpc, already mentioned in § |2.1| Believed to host a supermassive 
black hole of mass ~ 5 x lO^MofllEillit is well known for its low nuclear emissivity 
at all wavelengths. NGC 1399 exhibits weak FR I type radio activity (at a level of 
^ 10^^ erg/s; see reiW^ , but little evidence for significant non-thermal jet emission 
at higher energies. Based on an analysis of jet-induced cavities in the surrounding 
X-ray emitting medium on kpc-scalc, a jet kinetic power of Pbz ^ 10^^ erg/s has 
been inferred. ^^*^ If this jet power would be provided by a Blandford-Znajek-type 
(BZ) process, then ordered poloidal magnetic field strengths close to the black hole 
of the order of Bp ^ 300/a G (with a the spin parameter) are to be expected, cf. 
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eq. (281 



Estimates for the (ADAF-type) accretion rate based on Chandra observations, on 
the other hand, seem to suggest M ^ 0.1 A/Bondi ^ 2 x lO^^AfEdd (see refsEHMll) 
and would therefore imply (disk) magnetic field strengths (eq. 22) Ba ^ 10^ G, 



compatible with the one inferred from the observed jet power. Thus, if in addition 
the full electric potential would be available for particle acceleration (i.e., a ~ 1, 
h ^ Tg), electrons/protons might be able to reach maximum Lorentz factors ^ 10 



10 



i.e. protons may reach ultra-high energies of ~ 10""^^ eV, cf. eqs. ( |25| ) and (50 1. This 
would yield curvature emission in the ^ (0.1—1) TeV regime (eq. [51| ). An important 
consequence of this is that a system containing a rapidly spinning black hole (such 
that it can be an efficient UHECR accelerator) should also emit curvature TeV 
photons, provided (i) vacuum breakdown does not occur and (ii) the TeV photons 
are able to escapcl^^ In the case of NGC 1399 the later seems likely to be the case 
(cf. ref. 12231). 

In order to sustain the observed BZ jet power, a current / ~ (Phz/AZ^)^/^ ^ 
10^^ statampere ~ 3 x 10^^ A must flow through the magnetosphere. This would 
require an injection rate of dN^/dt ^ I /Q ^ 2 x 10^^ particles/s and imply a 
minimum charge density Ue ~ Nf./V ^ {dNf./dt)/{r'^c) ~ 10"^ particles/cm^ . If 
this cannot be provided by annihilation of MeV disk photons, an additional plasma 
source would be needed to estabfish a BZ-type jet. 

The accretion rate in NGC 1399 seems to be close to the critical value merit ^ 
3 X 10"'' where annihilation of MeV photons in an ADAF could lead to an injection 
of seed charges of density rie comparable to the Goldreich- Julian value ncj — 7 x 
10^^(S/300G) cm^3^ ^f^ ref-l^l If the accretion rate would be suflficiently high, i.e. 
TO > TOcrit ensuring rie > uqj, a substantial part of the electric potential may be 
screened (i.e., h < Vg likely) so that efficient gap- type particle acceleration becomes 
suppressed. The anticipated curvature VHE output would be of the order of (cf. 
eq. 52 ) Lc -- Pbz (h/rg)^. For {h/rgY ~ 0-1, this would be consistent with the VHE 
upper hmit L^{> 200 GeV) < 9.6 x lO'*'^ erg/s, imposed by H.E.S.S. observations 



of NGC 1399, see Fig. 12 ^^^ xhis would further reduce the available potential 
(eq. 26 ), limiting achievable proton energies to ^ 10^^ eV. If, on the other hand, 
the accretion rate would be sufficiently small (m < merit implying rig < tiqj)^ 
fully developed gaps {h ~ rg) may exist, but an additional plasma source (such 
as cascade formation in starved magnetospheric regions, ref.I^J) would be needed 
to ensure a force-free outflow. The gap would then emit VHE 7-rays with a total 
luminosity of order Lc ^ (rie/nGj) Pbz- For (n/ncj) < 0.1 this would again be 
consistent with the VHE upper limits from H.E.S.S. and Fermi observations. 

5.2.2. Quasar remnants as UHECR accelerators 

Provided screening is not effective (and the spin parameter high enough), accelera- 
tion of cosmic-rays to energies > 10^^ eV may become possible. This could perhaps 
be the case in inactive quasar remnants, harboring rapidly spinning, supermassive 
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Fig. 12. The spectral energy distribution (SED) of NGC 1399, compiled from non-simultaneous 
public data (VLA radio, ISO IR, HST, Chandra, Fermi and HESS). Above the optical regime, 
only upper limits are available. The thin black lin e rep resents a toy curvature spectrum assuming 
a fully developed gap. Figure from Pedaletti et al.EMl 



(> IO^Mq) black holes .l^^il 'j'jjg scenario envisaged is then one, where the black 
hole accelerator is not always operational in the normal, BZ-jet producing mode 
(associated with vacuum breakdown), but where (at least) sporadically the full gap 
size becomes available for the acceleration of a small number of protons. Primary 
target sources would thus be those that show little evidence for jets as well. While 
this may not work for NGC 1399 for reasons given above, Boldt & Ghosh^^'* ini- 
tially suggested a number of nearby massive dark objects that could potentially 
be efficient UHE (> 10^° eV) proton accelerators. However, when curvature losses 
are fully accounted for,l3S! most of these nearby sources turn out not to be capable 



of doing so, see eq. (50) (this seems also to apply to the revised list in ref! 



In fact, there is only a relatively small region of the parameter space (black hole 
mass, magnetic field strength) where one may expect proton acceleration to > 10^'^ 
to be marginally possible. For this parameter region, however, the total power of 
the accompanying VHE radiation would appear to exceed the total power emitted 
in UHE cosmic rays by a factor of a hundred or more, making such hypothetical 
sources no longer quiescent in the TeV regime.'SSl Given current evidence, it seems 
thus rather unlikely that (nearby) massive quasar remnants are promising candi- 
dates for efhcient (gap-type) acceleration of protons to energies beyond 5 x 10^^ 
eV, even if energy losses due to photo-pion production, cf. ref.'2251, would be ne- 
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glected. As the constraint imposed by curvature losses is only mildly dependent on 
the charge number, heavier nuclei such as, e.g. iron Fe, might in contrast be able to 
reach energies ~ 5 x 10^° eV in massive remnants. The propagation of such nuclei 
will, however, be significantly influenced by the interactions with the intergalactic 
radiation fields through photo-disintegration, with the particle being stripped from 
all its constituents nucleons after distances of a few Mpc onlyEMUZl 

6. Conclusions 

Non-thermal processes occurring in the vicinity of supermassive black holes may 
allow new insights into the behavior of matter and energy under extreme condi- 
tions. The detection of, e.g., rapid variability on timescales of ~ r^/c and unusual 
high-energy emission signatures can provide an important diagnostic tool towards 
this end. Based on this, VHE observations during the last years have demonstrated 
the capacities of the new generation of VHE gamma-ray instruments (e.g., MAGIC 
II, HESS II, VERITAS and FERMI) to substantially improve our knowledge of the 
activity center in nearby, underluminous AGNs and to trigger new theoretical devel- 
opments that could serve as benchmark for future VHE arrays like CTA. Detailed 
theoretical modelling of non-thermal particle acceleration and radiative numerical 
simulations may soon become indispensable in order to fully take advantage of the 
improved observational capacities. 
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